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ABSTRACT 


If the current growth rate in the use of fossil fuels continues 
at 4.3% per year, then the CO^ concentration in the atmosphere can 
be expected to double by about 2035 provided the current partition 

C °2 between the atmosphere, biosphere and oceans is maintained as 
is the current mix of fuels. Slower rates of anticipated growth of 

energy use lead to a doubling of the carbon content of the atmos¬ 
phere sometime in the period 2040 to 2060. 

This report addresses the questions of the sources of atmospheric 
carbon dioxide, considers distribution of the present carbon dioxide 
among the atmospheric, oceanic and biospheric reservoir and assesses 
the impact on climate as reflected by the average ground temperature at 
each latitude of signfleant increases in atmospheric carbon dioxide. 

The oceans and the biosphere, including the soils, store very large 
amounts of carbon compared with that in the atmosphere. Small changes 
in these large reservoirs can have a major effect on the atmosphere. 

Neither the ocean-atmosphere nor the atmosphere-biosphere interfaces 
are well understood. In regard to the former we propose a new model 
for the mixing of carbon dioxide in the oceans. The proposed model 

explicitly takes into account the flow of colder/and or saltier water 
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to great depths. Such models provide a means for the evaluation of the 
oceanic carbon dioxide mixing without resort to the use of box model 
which with adjusted effective diffusion parameters. The atmosphere- 
biosphere interface is significantly more complex. It is clear, 
however, that data on several key questions are needed: (1) clear and 
detailed information on the rapid turnover time of life forms in the 
oceans (2) information on the possibility that biologically derived 
organic debris in shallow seas could be buried and thus effectively 
removed from the oceans (3) the reservoir of carbon in the soils is 
thought to be about five times the amount in the atmosphere but the 
fraction of that carbon which can communicate with the atmosphere is 
not known. 

Increasing CO^ in the atmosphere perturbs climate by altering 
the radiative properties of the atmosphere. The resulting climate 
change has been calculated through highly parameterized analytic models 
of the ocean-atmosphere system or by large computer models. We have 
adopted the first approach starting with an atmosphere locally in radia¬ 
tive equilibrium and then allowing meridional transport of heat. We 
have constructed two models for the case of radiative equilibrium 
treating the atmosphere as gray and dividing the infrared emission 
region into nine bands. The gray atmosphere model predicts an 
increase of average surface temperature of 2.8°K for a doubling of 
CO^, a result about a degree less than the nine band model. 

In the model the principal effect of increasing CO^ is to enhance 
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the absorption by weak CO^ bands in 8-12 micron region. Trace gases, 
CH^, N^O, NH^, freons and hydrocarbons can also block off the at¬ 
mospheric window. 


An analytic model of the atmosphere has been constructed (JASON 
Climate Model). Calculation with this zonally averaged model shows 
an increase of average surface temperature of 2.4° for a doubling of 
C0 2 . The equatorial temperature increases by 0.7°K while the poles 
warm up by 10 to 12°K. 


The JASON climate model suffers from a number of fundamental 
weaknesses. The role of clouds in determining the albedo is not ade¬ 
quately taken into account nor are the asymmetries between the northern 
and southern hemisphere. We expect, however, that models intermediate 
between the large GCMs and the primitive analytic models can yield 
insights into the nature of climate change. 


The warming of climate will not necessarily lead to improved 
living conditions everywhere. Changes in sea level, in agricultural 
productivity and in water availability can be anticipated, but the 


dimensions of their economic, political or social consequences can not. 
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illuminated many of the problems of modeling the atmosphere for long¬ 
term climate changes. Their perceptive presentations convinced us to * 

go back to basic principles and try to understand radiative equilibrium 
and energy budgets. We have tried. 
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As physicists, we welcomed the opportunity to provide the 
rationale that biology was not in any way of significance in under¬ 
standing effects of carbon dioxide. ^ 
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convinced us that the biosphere plays a critical role in the carbon 
dioxide problem. We are now attempting to understand the carbon 
dioxide consequences of life in the oceans with its very small standing 
crop and very rapid turnover, and life on land where opposite conditions 
hold. 
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THE IMPACTS OF INCREASING CARBON DIOXIDE: AN OVERVIEW 


Events in the early 1970's illustrate the influence that small 
shifts in climate can have on world events. In 1972 the Soviet Union 
lost a significant portion of its winter wheat crop when the snows 
failed to provide adequate cover, and a sharp freeze destroyed the 
exposed seedlings. Later, the winter harvest was delayed as the summer 
moisture normally carried by the westerlies did not arrive in the Ukraine 
and in the late fall, the ripened crop was blanketed by snow. The 
Soviet Union rapidly entered the world grain market with major world¬ 
wide economic repercussions. 

In the same time interval, the six West African countries south 
in the Sahel—Mauretania, Senegal, Mali, Upper Volta, Niger and Chad, 
were driven to the edge of economic and political ruin by a drought 
that began in 1968. The Sahelian drought and the Soviet grain failure 
received world-wide attention. The Ecuadorean rice failure of 1974, the 
Chinese droughts and floods of the same year and the 1973 droughts in 
Central America attracted less attention. However, all of these 
examples illustrate the fragility of the world's crop producing 
capacity, particularly in those marginal areas where small alterations 
in temperature and precipitation can bring about major changes in 
total productivity. 
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In 1977 and 1978, climatic conditions for agriculture largely 
returned to those of the 1950s and 1960s with agriculturally favorable 
weather prevailing in the Indian subcontinent, China, the Soviet Union, 
and the United States, Despite the improved world food situation, the 
lesson of the early 1970 f s remains, the heightened demands of an in¬ 
creasing population for food with a lowered availability of arable land 
and water make the world food supply vulnerable to climatic changes. 

The climatic fluctuations of the past few decades are minor compared 
to those contained in the historical and geologic records. A warming 
comparable to the "altithermal" period of 4000 to 8000 years ago would 
bring about not only major shifts in agriculture, but would have sig¬ 
nificant demographic effects not only from the direct climatic changes, 
but also from shifts in sea level resulting from the melting of glacial 
ice. 


The burning of carbon-based fuels employing current technology 
releases carbon dioxide into the atmosphere. Carbon dioxide is a 
relatively minor constituent in the atmosphere, having a concentration 
at present of 332 parts per million (ppm). Despite its small con¬ 
centration, carbon dioxide exerts a major influence on the thermal 
structure of the atmosphere since the carbon dioxide molecule absorbs 
infrared radiation emitted by the earth’s surface that otherwise would 
escape into space. Because of carbon dioxide’s radiative properties, 
changes in the concentration of carbon dioxide changes the atmosphere's 
thermal properties, and in this fashion changes in carbon dioxide 
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concentration can produce changes in climate. While the general theory 
of the influence of carbon dioxide on climate is widely accepted by 
atmospheric scientists, there are disputes as to the magnitude and 
character of the climate fluctuations resulting from anticipated future 
burning of carbon based fuels. 

Plants in photosynthesis fix carbon from carbon dioxide in the 
atmosphere. Changes in carbon dioxide concentration can be expected 
to alter biological productivity, but in ways that are not understood. 
The acidity of the buffered oceans depends in part on the atmospheric 
carbon dioxide concentrations and changes in this concentration could 
influence marine life. These examples emphasize that while climatic 
variations are likely to be the most important consequences of changes 
in carbon dioxide, other changes whose importance is presently not 
known, can be anticipated. 


A. Climate and Carbon Dioxide in Historical Perspective 

Tyndale^ in 1861 first suggested that slight changes in atmospheric 

composition could bring about climatic variation. Tyndale’s work was 

based on early experiments measuring the absorption of light as a 

2 

function of gas density. Arrhenius first calculated the influence of 
changes in carbon dioxide concentration on the surface temperature noting 
that enhanced absorption in the infrared should yield a higher surface 
temperature. For a three-fold increase in the atmospheric concentration 
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of carbon dioxide, Arrhenius in 1896 calculated a surface temperature 

rise of 9°C—a remarkable result considering the paucity of data. 1 

3 4 

Today’s estimate is 4 to 6°C. Arrheniys and Chamberlain 
speculated that the large variation in the earth’s climate, and in 

particular the glacial epochs, could be due to changing carbon dioxide 4 

concentrations of the atmosphere resulting from geologic processes. 

The work of Arrhenius and Chamberlain greatly influenced the thinking 

about geologic changes in climate during the first half of the 20th 4 

Century. The critical role of the oceans in the global distribution 

of carbon dioxide was first demonstrated and explored by Tolman^ in 

1899. 4 

Calendar in 1938 first recognized that man, through the burning 
of carbon-based fuels, was changing the composition of the atmosphere. 4 

The geophysical significance of the rise in atmospheric carbon dioxide 
was dramatically emphasized by Revelle and Suess^; "Human beings 

are now carrying out a large-scale geophysical experiment of a kind 4 

that could not have happened in the past nor be repeated in the future. 

Within a few centuries, we are returning to the atmospheres and oceans 
the concentrated organic carbon stored in the sedimentary rocks over to 

hundreds of millions of years. This experiment, if adequately docu¬ 
mented, may yield a far-reaching insight into processes determining 

weather and climate." I 

Revelle was instrumental in incorporating into the program of 
the International Geophysical Year, accurate and regular measurements t 
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has led a number of 


8 9 

of the concentration of carbon dioxide. Keeling * 
investigators in maintaining, over the past 20 years, continuous monitoring 
programs of atmospheric carbon dioxide at Mauna Loa, Hawaii, and the 
South Pole. 

The impact of rising carbon dioxide levels on man were alluded to 
by Revelle and Suess and explicitly discussed in a report of the 
President's Science Advisory Committee^ in 1965. The Arab Oil 
Embargo, President Nixon's Project Independence and the increased 
emphasis on the use of coal as a major energy resource has focused re¬ 
newed attention on the possible consequences of increased levels of 
carbon dioxide in the atmosphere. 


B. Measured Changes in Atmospheric Carbon Dioxide 

Figure 1-1 shows the record of a secular increase of atmospheric 

carbon dioxide obtained from a nearly continuous monitoring at Mauna 

Loa, a site at 3,400m altitude and well within the trade wind belt. 

The current content carbon as carbon dioxide in the atmosphere is 
9 

about 702 x 10 metric tons (Gtons). Between about 87 and 140 Gtons 
of carbon has been added to the atmosphere since 1850 (See Table 1-1). 
The uncertainty arises from the unknown level of carbon dioxide in 
the pre-industrial atmosphere (See, for example, Callendar^). 
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Figure l-l 

Atmospheric Carbon Dioxide Concentration 
at Mauna Loa Observatory 

(after Keeling, et el 8, 9) 




TABLE 1-1 


CHANGING C0 2 CONTENT 
OF ATMOSPHERE 



CONCENTRATION 



IN PPM BY 

MASS OF CARBON 

DATE 

VOLUME 

IN Gtons) 

1850 

265 - 290 

562 - 615 

1958 

315.A 

669 

1975 

330.9 

702 
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A similar secular increase has been observed at the South Pole 
with the values lagging Mauna Loa by a few ppm in any given year. This 
lag is consistent with the interpretation that more than 90% of the 
increased carbon dioxide is introduced in the industralized Northern 
Hemisphere, and that there is a time lag in the mixing of the atmosphere 
between the hemispheres. 

Shorter records from Point Barrow, Alaska^ and from Scandinavin 
12 

airplane flights show similar secular trends. Superimposed on the 

secular trend is a regular annual variation of approximately 6 ppm 

at Mauna Loa, 1.6 ppm at the South Pole, and as high as 15 ppm at 

Point Barrow. In the Northern Hemisphere, the maximum in the annual 

variation is in spring and the minimum in late summer. The seasonal 

variation is usually attributed to seasonal changes in net photo- 

13 

synthesis and respiration of the biosphere though variations in sea 
surface temperature may contribute since the dissolved carbon dioxide 
content of the oceans is temperature dependent. 

Even with the limited number of observation stations, it is 
clear that the carbon dioxide content of the atmosphere has been rising 
and that the rate of increase has accelerated over the 20 years of 
accurate measurements. In the late 1950 , s, carbon in the form of 
carbon dioxide was being added at a rate of about 1 Gton per year, 
while today the rate is approximately two Gtons per year. 
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c. 


Fuel Sources of Carbon 


Calender** first suggested that the burning of carbon-based fuels 
provided the major source for the increasing carbon dioxide content of 
the atmosphere. The burning of calcium carbonate to make quick lime 
in cement manufacture also adds carbon dioxide, but this contribution 
of carbon dioxide returns to cement once the cement and mortar are used. 
Table 1-2 lists the carbon added to the atmosphere by the world-wide 
burning of carbon-based fuels for various years. Over the last 30 years, 
the fuel generated carbon has been increasing at an almost constant rate 
of 4.3% per year. 

The cumulative contribution by fuel burning of carbon to the short¬ 
term carbon cycle for the period 1800 to 1978 is about 154.5 Gtons or 
about one-quarter of the pre-Industrial atmospheric carbon content. 
However, as will be discussed, only a fraction of the carbon introduced 
into the atmosphere has remained airborne. Of the 154.5 Gtons, 27% or 
42.5 Gtons have been added in the period 1970-1978. 

The amount of carbon added to the atmosphere depends on the mix of 
carbon-based fuels. Table 1-3 lists values for carbon added per hundred 
quads (1 quad - 10'*'"* Btu) of thermal energy generated. The values are 
estimated averages since the heat content of the fuels vary, particularly 
for coal. The values for synthetic gas and oil derived from coal are 
even more uncertain because of the differences among possible gasifica¬ 
tion and liquefaction processes. The important point is that coal derived 
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TABLE 1-2 


ESTIMATED CARBON ADDED TO THE ATMOSPHERE 
BY THE BURNING OF FUELS IN GTONS 14 


YEAR 

CARBON ADDED 

1950 

1.63 

1960 

2.61 

1970 

3.96 

1975 

4.87 

1978* 

5.54 

*Estlmated 



TABLE 1-3 

CARBON PRODUCED PER HUNDRED QUADS OF 
THERMAL ENERGY GENERATED 

CARBON 

FUEL IN GTONS 

Oil 
Gas 
Coal 

Synthetic Gas 
Oil Derived 
From Coal 


2.0 
1.45 
2.5 
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TABLE 1-4 


1975 WORLD PRODUCTION OF CARBON 
BY BURNING FOSSIL FUELS 


FUEL 

QUADS 

PERCENT 

GTONS 

OF CARBON 

Oil 

112 

47% 

2.24 

Gas 

49 

20% 

0.71 

Gas (flared) 

10 

4% 

0.14 

Coal 

70 

29% 

1.78 

TOTAL 

241 

100% 

4.87 

4.87 Glgatons 

- 0.69% of 

total atmospheric carbon 
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synthetic fuels will release more carbon dioxide per unit of delivered 
thermal energy than will coal. The relative mix of carbon-based fuels 
and their contribution in 1975 to the atmospheric carbon load is shown 
in Table 1-4. 

I 

Projection of future carbon loading of the atmosphere is clearly 
a hazardous venture depending on anticipated rates of energy usage in 

various parts of the world as well as on estimates of the mix of fuels • 

employed by various countries. For example, the Soviet Union plans to 
increase the fraction of their total energy produced by natural gas 

from the present 22% to 40% in 1985"^. The switch from coal and oil • 

to gas will lower the amount of carbon introduced into the atmosphere 

per quad of energy used. Similarly, large scale introduction of nuclear 

or solar energy systems will further reduce the rate at which carbon • 

is introduced. 

If the use of carbon-based fuels continues to grow at a yearly 4 

rate of 4.3%, then 14.4 Gtons of carbon would be released as carbon 
dioxide in the year 2000, 34 Gton in 2025. For a continuing 4.3% rate 

of increase the cumulative total of carbon released by 2000 would be 4 

about 360 Gtons, by 2025, 840 Gtons. These estimates would probably 
be upper limits since energy growth rate will slow in the industrial 

world and the use of non-carbon fuels will increase. However, the 4 

extensive use of coal in the next 20 years could raise the yearly 
contribution somewhat (See Table 5). 

4 
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TABLE 1-5 


WORLD-WIDE 

CARBON RELEASE 

IN THE YEAR 

2000 

WORLD CONSUMPTION 

FUEL MIX 


OF CARBON BASED 

IN PERCENT 

CARBON RELEASE 

ENERGY IN QUADS 

OIL GAS 

COAL 

IN GTONS 


50 20 

30 

12.2 

600 

30 10 

60 

13.5 


30 50 

20 

10.9 



50 

20 

30 

14.0 

700 

30 

10 

60 

15.8 


30 

50 

20 

12.7 


TABLE 1-6 

WORLD-WIDE CARBON RELEASE IN THE YEAR 2025 


WORLD CONSUMPTION 


FUEL MIX 



OF CARBON BASED 

IN PERCENT 

CARBON RELEASE 

ENERGY IN QUADS 

OIL 

GAS 

COAL 

IN GTONS 

770 

20 

0 

80 

18.4 

10 

60 

30 

14.0 

1150 

20 

0 

80 

27.5 

10 

60 

30 

20.9 
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Table 1-5 presents Illustrative values of carbon release based on 


« 

a low (600 quads) and high (700 quads) yearly of energy uses in the • 

year 2000. Depending on the assumed carbon fuel mix the values range 
from twice to three times the current release rate. Similar estimates 
are shown in Table 1-6 for the year 2025 where the annual energy use rate It 

of 770 quad corresponds to a 1% yearly growth from a base of 600 quads 
in 2000, the 1150 quad figure corresponds to a 2% growth rate in use of 
carbon-based fuels from a base of 700 quads in the year 2000. The car- t 

bon release 45 years in the future could thus be three to six times the 
current rate and even greater if synthetic fuels are used in a signifi¬ 
cant way. « 

D. The Carbon Mystery f 


A comparison of the observed increase of atmospheric carbon 
(Figure 1-1) with the release of carbon from burning fuels (Table 1-4) |i 

shows that only about half of the fuel carbon remains airborne and thus 
capable of influencing climate. The usual assumption is that the oceans 
absorb the majority of the missing carbon though a portion may be taken t 

up by the biosphere^. More recently, the possibility has been advanced 
that current land clearing with the subsequent oxidation of plant 
- material and the removal of carbon using plants may be adding to the 14 

carbon content of the atmosphere in a significant way. An understanding 
of the distribution of carbon between the oceans, biosphere, sediments 
and atmosphere is of critical importance to the forecasting of future t 

carbon levels in the atmosphere. 
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TABLE 1-7 

APPROXIMATE MAGNITUDES OF THE 
RESERVOIRS OF CARBON IN GTONS 


• Sediments 

- Oxidized 5 X 10 7 

- Reduced 2 X 10 7 

• Igneous Rocks 3 x 10^ 

4 

• Oceans (deep) 4 X 10 

2 

• Oceans (mixed layer) 6 X 10 

• Fuels 7 X 10 3 

• Soils 3 X 10 3 

• Biosphere (living) 8 X 10^ 

2 

• Atmosphere 7 X 10 
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The atmosphere is in fact the smallest of the major reservoirs 
of carbon as is illustrated in Table 1-7. Further, the fluxes between 
the reservoirs are large. For example, the biosphere puts into the 
atmosphere approximately 160 Gton of carbon a year through respiration 
and decomposition of terrestrial and oceanic biota. The atmosphere 
returns an equal quantity through photosynthesis if equilibrium 
conditions are maintained. Small changes in the large reservoirs or 
in the rates of exchange between the reservoirs could bring about 
significant changes in atmospheric composition that would overwhelm the 
effects of burning carbon based fuels. Indeed, the older climate 
change hypothesis put forward by Arrhenius and Chamberlain were based 
on changes in the carbon content of the various reservoirs on geologic 
time scale. 

The uncertainty in how much of the man-derived carbon goes 
rapidly into the oceans derives from the poorly understood mixing 
processes within the ocean. Numerous models of ocean-wide diffusive 
and convective exchange processes between adjoining vertical layers 
have been constructed to study the movement of carbon dioxide. The 
models and their supporting data on the movement of radioactive iso¬ 
topes are not accurate enough to determine what fraction of the carbon 
released through fuel burning enters the ocean or whether additional 
carbon derived from land clearing and deforestation is also sequestered 






« 
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In the oceans. In this report, Munk and Zachariasen (§ II-B) question 
the application of convection-diffusion models arguing that the movement 
of carbon dioxide in the oceans is influenced by processes at the ocean 
boundaries which alter the temperature and salinity (hence the density) 
and which affect how the newly formed water masses spread into the oceans. 

The question as to whether the biosphere is a source or sink for 

carbon is similarly uncertain. The rates of land clearing, the resulting 

changes in biological productivity, the removal of organic material by 

sedimentation in the shallow oceans, the rates of exchange of carbon 

in soils with the atmosphere are all poorly known. Recent estimates 

suggest that the biosphere has been a moderate source of carbon in 

18 

recent years, 1 to 3 Gtons/year , but much higher figures have been 

proposed as well as suggestions that the biosphere is a sink for carbon. 

Some oceanographers would argue that the ocean data are inconsistent 

with the hypothesis that the oceans take up half of the carbon supplied 

by fuel burning (2.5 Gtons/year) together with the 1 to 3 Gtons/year 

19 

derived from changes in the biosphere. 


E. Future Carbon Levels in the Atmosphere 

Given uncertainties in future carbon fuel use, fuel mix, biospheric 
contribution and ocean-atmosphere interaction, it is hazardous to pre¬ 
dict future levels of carbon dioxide in the atmosphere. If the present 
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rate of increase (4.3% per year) of use of carbon based fuels continues 
indefinitely into the future, then the carbon dioxide content of the 
atmosphere would double by the year 2035. Table 1-8 illustrates the 
dependence of the date of doubling of atmospheric carbon dioxide on 
the rate of energy use assuming the current mix of fuels, and that half 
of the carbon remains in the atmosphere. At conservative rates of 
increase of carbon fuel use (one to two percent per year), assuming 
that half of the fuel carbon remains airborne, the carbon dioxide 
levels would be expected to double current levels (332 ppm) in 60 to 
80 years. The time would be reduced if future agricultural pressures 
lower the standing crop, or if coal use, particularly synthetic, in¬ 
creases at a greater rate than use of oil or gas. 

Since large reserves of carbon-based fuels are known to exist, and 
much larger resources are suspected (see Tables 1-9 and 1-10), the burning 
of the fuels could increase the carbon content of the atmosphere five 
or ten fold over a time span of centuries if carbon-based fuels continue 
to be significant energy source. 

F. Climate Change From Changing Carbon Dioxide 

Attempts to calculate the effects of carbon dioxide on climate 
have focused on determining the change in the average surface temperature 
of the earth. Climate is much too complicated to be described by a 
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TABLE 1-8 


APPROXIMATE DATES FOR DOUBLING THE 
CARBON DIOXIDE CONTENT OF THE ATMOSPHERE 


ASSUMED RATE OF 
INCREASE IN USE OF 
CARBON FUELS 


Constant 4.3% 
Constant 3% 

4.3% 1978 - 2000 

3% 2000- 
2 % 2000 - 
1 % 2000 - 


DATE OF DOUBLING 
OF CARBON DIOXIDE 
CONTENT OF ATMOSPHERE 

2035 

2050 


2040 

2050 

2060 
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TABLE 1-9 


PROJECTED WORLD CARBON PRODUCTION 
BY USE OF RECOVERABLE RESERVES 20 

RECOVERABLE 



RESERVES 

CARBON RELEASF 

FUEL 

(QUADS) 

(GTONS) 

Oil 

3,200 - 3,700 

64 - 74 

Gas 

1,700 - 2,600 

24 - 37 

Coal 

15,000 

381 

TOTAL 

19,900 - 21,300 

469 - 492 


TABLE 1-10 


PROJECTED WORLD CARBON PRODUCTION 
BY USE OF RECOVERABLE RESOURCES 20 


FUEL 

RECOVERABLE 

RESOURCES 

(Quads) 

CARBON RELEASE 
(GTONS) 

Oil 

8,700 - 11,000 

174 - 220 

Gas 

9,100 - 9,800 

131 - 141 

Coal 

86,000 - 160,000 

2,184 - 4,064 

Sub-Total 

103,800 - 180,800 

2,489 - 4,425 


Oil shales and 
Tar sand 

Non-conventional 
Gas (U.S. only) 

TOTAL 


14,000 - 100,000 
20,000 - 60,000 

137,800 - 340,800 


280 - 2,000 
288 - 864 

3,057 - 7,289 


20 


Digitized by 



single parameter. Amounts and time distribution of precipitation, onset 
of freezing conditions, strength and patterns of wind are all parameters 
critical for understanding the impact of climatic change on man. However, 
present knowledge is such that surface temperature is universally used 
as a surrogate for climate in calculation on climatic change. 

Three principal classes of assumptions have been employed to 
calculate the effects of carbon dioxide on surface temperature. 

1) Assume the atmosphere is everywhere in radiative equilibrium 
with outgoing infrared radiation balancing Incoming solar radiation at 
each point on the earth*s surface. (Zeroth order model) 

2) Assume that in addition to radiative processes, heat is 
transported meridionally from equator to pole in the oceans and 
atmosphere by eddy-diffusion. (First order model) 

3) Assume that the basic hydrodynamical equations describing 
atmospheric motion can be solved numerically. (General circulation 
models, GCM) 


G. Zeroth Order Calculations 


The earth*s thermal radiation is mainly confined to the 5 to 30 pm 

region. The water molecule is a strong absorber over much of this 
region, except for the important 8 to 18 pm interval where the earth's 
infrared emission peaks (for a surface temperature of 290°K, the peak 
is at a wavelength of 10 pm). Water is a principal contributor to the 
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greenhouse effect, even though it is not everywhere optically thick. 

Climate calculation usually assumes that the water vapor concentration 
remains constant even though carbon dioxide and temperature vary. 

Current carbon dioxide abundances are enough to make the 15 pm 
region optically thick so that changes in carbon dioxide concentration 
will alter the greenhouse heating primarily by carbon dioxide absorption 
in the 10 pm band. The situation is further complicated in that a 
number of minor constituents such as nitrous oxide, methane, and a wide 
array of hydrocarbons are potentially important In enhancing the effect 
of carbon dioxide in closing off the 8 to 15 pm region. For example, 
a doubling of the methane concentration, currently 1.6 ppm would pro¬ 
duce a temperature increase in excess of 1°C. 

Chamberlain (§ ill-A) has calculated on the basis of radiative theory, 
that if the carbon dioxide were doubled, then the average surface 
temperature would Increase by 2.8°C using the assumption of a gray homo¬ 
geneous atmosphere. Vesecky (§ III-B), using a nine-band model obtains 
a temperature increase of 3.8°C. As carbon dioxide is increased 
further, the temperature will continue to rise, although the 15 pm bands 
become more saturated, the 10 pm bands are effectively closing the 8-12 pm 
window. Chamberlain's and Vesecky's recent calculations are in general 
agreement with earlier radiative calculations. 
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H. First Order Models 


In the first order models, the surface temperature is permitted 
to vary with latitude. Heat transport is modeled as a eddy-diffusive 
process and the oceans and atmosphere are considered to be a coupled 
system. The resulting equations are non-linear since the albedo is 
assumed to be a function of temperature. In the simplest case, the 
albedo is high where the temperature is low because of ice or snow 
cover. Carbon dioxide will affect the radiative components in the 
equation and thus influence both the average temperature and the 
temperature distribution with latitude. 

Abarbanel (§ IV-B), using the First Order Model approach, finds that 
doubling the carbon dioxide content from 332 ppm to 664 ppm results in 
a temperature increase of 2.4°C. Interestingly, he finds that the 
temperature rise is much more marked in the polar regions with a 
temperature increase of 10°C at high latitudes, while in the equatorial 
regions it is less than the average. 


I. Global Circulation Models 


20 

Manabe and Wetherald have carried out the most detailed computing 
machine calculation using numerical simulations of the atmosphere but 
of necessity approximating very roughly cloud formation, air-sea inter¬ 
action, and albedo change. They found that a doubling of the carbon 

23 


Digitized by v^ooQle 


dioxide content raises the average temperature 2.9° and the calculations 
indicate a warming of 8 to 20°C in the high latitude regions with 
smaller than average increases in the equatorial regions. 

A wide variety of models lead to similar numerical results with 
respect to the effect of doubling the carbon dioxide content, as has 
been noted by Schneider^. This observation is suggestive that the 
real control on surface temperature is due to the effects on the 
radiative energy balance by carbon dioxide and that the detailed 
dynamics of the atmosphere may not be crucial in determining average 
changes. A number of calculations are also in general agreement with 
the average warming being amplified at high latitudes. 


J. Impacts of Climate Change 

If civilization continues its heavy reliance on carbon based fuels, 
and if there are no major shifts in the current response of the oceans 
and biosphere to changing carbon dioxide content, then we should expect 
during the middle of the 21st Century a warming of 2 to 3°C accen¬ 
tuated by a factor of three or four at high polar regions. These 
average temperature changes are greater by at least a factor of two 
than those observed during abnormal historical weather conditions that 
prevailed for example during the M Little Ice Age' 1 of the 16th and 17th 
Centuries. The changes, however, are smaller by about a factor of two 
than those experienced during the past million years in the glacial 
and interglacial periods. 
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The anticipated man-made changes by the year 2000 or even 2025 are 
likely to be small compared with historical variations though carbon 
dioxide, together with minor constituent pollutants, could either 
amplify a natural warming trend or dampen a natural cooling episide. 

At the present state of knowledge, we cannot predict whether the natural 
variation will be in the direction of warming or cooling. 

Viewed in the longer perspective and noting past historical events, 
changes in climate would not lead to improved living conditions every¬ 
where. Any change would produce stress and possibly disaster in some 
parts of the world since so many aspects of society have adapted, with 
very large investments in their infra-structures, to the climate of the 
past few decades. 

Perhaps the most ominous feature of an average warming due to 

increased carbon dioxide, is the amplification of the warming in the 

polar regions. Generally, it is assumed that the polar ice sheets 

would respond to a warming trend over a time scale of thousands of years 

since the melting would proceed slowly at the surface exposed to the 

22 

atmosphere. Mercer , however, has pointed out that the characteristics 
of the West Antarctic ice sheet are such that melting could proceed 
rapidly. The ice sheet is grounded below sea level with present summer 
surface temperatures of -4 to -5°C. The grounded portion is fringed 
by ice shelves which serve to buttress the grounded part of the ice 
sheet. The shelves are vulnerable to both oceanic and atmospheric 
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warming. An increase in summer temperature of 0 to 1°C could lead to 
rapid melting of the shelves and disintegration of the main ice sheet. 
The time scale for such a process is not known, but observed surges 
indicate that the West Antarctica ice sheet could discharge one-third 
to one-half of its volume in 100 years. The complete melting of the 
ice sheet would raise world-wide sea levels by 5m with consequent major 
disruptions of the world’s coastal region, A warming of 5°C in 
Antarctica could take place in 50 to 75 years, depending on rates and 
kinds of carbon fuels burned. 




m 


% 




A warmer climate will impact agriculture primarily through 
lengthening the growing season in higher latitudes and increasing 
summer temperatures provided precipitation patterns are not greatly 
altered. Higher summer temperatures will not necessarily increase 
productivity of all crops, since the increase of respiration could 
exceed the increases in photosynthesis. In high latitude regions, 
a northward movement of a viable growing season does not imply greater 
productivity, since the soils on formerly glaciated terrain are 
generally poor. 

Some historical evidence suggest that warming could result in 
prolonged mid-continental drought' 1 ' 3 . n Dust Bowl" conditions could 
threaten range lands and semi-arid agriculture over large areas of 
North America, Asia and Africa. 
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Warming could also affect water availability both through altering, 
in unknown ways, precipitation patterns and enhancing evaporation from 
soils and plants. The latter effect, in decreasing soil moisture, 
reduces the amount of water available for groundwater reserves. 

Warmer climates will increase evaporation from lakes and reservoirs. 
Thus, one possibility is that climatic warming could alter the water 
balance over large regions. 

The net energy impact of a warmer climate is not known. Demands 
on heating fuels will decrease, but these would be accompanied by 
increased use of air-conditioners. The net balance would be influenced 
by the demographic shift. Induced by climatic change. 

Shifts in coastlines, in agricultural patterns and water avall*- 
ability are sure to have major demographic effects. The social and 
political consequences of shifts in patterns of population are surely 
uncertain, but could be highly disruptive. 


K. Conclusions 


Despite the many uncertainties about the nature of the carbon 
cycle and of the effect on carbon dioxide on climate, it seems highly 
probable that continued increased in the world-wide use of carbon 
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based fuels at a rate of one to three percent per year, will lead to 
climatic changes in the second half of the 21st century. The changes 
will be amplified in high latitude regions and in these regions man¬ 
made changes may be larger than natural variations in the first half 
of the 21st Century. The impact of the projected changes on man is 
unknown but unlikely to be wholly favorable. 




« 


Means of ameliorating the effects of carbon dioxide increases 
are available, but application of these are unlikely to alter the 
above conclusion except for lengthening the time scale over which 
expected climatic changes would take place. The increased use of 
natural gas, if available, would lengthen the time available for a 
shift to non-carbon based fuels. Clearly a significant shift to a 
nuclear or solar energy economy would postpone carbon-induced climate 
shifts. Increasing the standing crop through massive reforestation 
could provide temporary storage for carbon, but land and water avail¬ 
ability limit this measure. In principle, carbon dioxide in stack 
gas could be frozen into dry ice for deep ocean disposal or the stark 
gases could be used to accelerate the growth of certain plants in 
giant greenhouses. However, elementary considerations suggest that 
neither of the latter suggestions are economical in competition with 
alternative, nuclear or solar, energy systems. 
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The potential problems posed by a warming climate appear 
sufficiently serious to justify a comprehensive research effort de¬ 
signed to reduce the many uncertainties discussed above. Even though 
anticipated alteration would appear to be at least 50 years in the 
future, the world-wide nature of the changes and their possible adverse 
effects on man-kind warrant the continued attention of policy-makers 
to the carbon dioxide climate question. 
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II 


THE OCEANS AS A SINK FOR CARBON DIOXIDE 


The atmosphere interacts with the oceans through diffusive pro¬ 
cesses at the air-sea interface. Carbon dioxide mixes with the sea 
water and the carbon first resides as a dissolved gas or as a component 
of the carbonate or bicarbonate ion. The carbon can then be taken up 
by living matter to become organic carbon. The details of the many 
pathways by which carbon can move through the ocean is poorly under¬ 
stood since chemical, biological and physical processes are all Involved. 

Models of the ocean have been constructed in which the ocean is 

divided into boxes. The early models divided the ocean into two layers, 

an upper layer called the mixed layer, in which turbulent diffusive 

processes even out the carbon content which is in approximate equili- 

1 2 

brlum with the atmospheric carbon * . The deep layer exchanges carbon 
with the mixed layer, but the two layers are not at any one time in 
chemical equilibrium. 

14 

Studies of the Suess effect, radioactive C dilution by the 

14 

burning of carbon based fuels, bomb injected C and other radioactive 

3 4 5 6 

isotopes have led to box models of the oceans of increasing complexity. 9 9 9 
These models are limited in the sense that data on the distribution of 
the isotopes must be used to determine the empirical transport 
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coefficients so that all one can hope for is a consistency between the 
observations and the model results. The box models by themselves cannot 
illuminate the physical, chemical and biological mechanisms by which 
carbon travels from the atmosphere through the oceans to sediments 
or back to the atmosphere. Because of the limitations of the box 
model, we explore in a preliminary way alternatives that use oceano¬ 
graphic data on the mass flow within the oceans. 


A. Cold Ocean Sinks 


The air-sea interface at which CO^ is enchanged between the at¬ 
mosphere and the oceans is adjoined to a well-mixed sea surface layer 
whose average thickness is conventionally taken as about 75 m. The 
total carbon in the layer (about 580 Gtons) is very comparable to 
that present as CO^ in the atmosphere (702 Gtons); but only 1% of 
the mixed layer carbon is in the form of dissolved CO^ - 90% is in 

the carbonation CO^ and almost all of the rest is HCO^. The rela- 

2 

tively small mixed layer volume thus contains about 10 times the 
carbon density of the adjacent atmosphere by storing almost all the 
excess in carbonate and bicarbonate ions. But the layer is, at present, 
about 85% saturated and it is, consequently, difficult for it to 
absorb much more: almost all additional absorbed CO^ must be accompanied 
by the disappearance of a carbonate ion effectively through the reac¬ 
tion CO^ + CO 2 + H^O - ZHCO^. Because of the large amount of CO^ 
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already absorbed, there is little CC> 3 still available and further 
increases in atmospheric would be accompanied by very much smaller 
increases in the carbon content of a mixed ocean layer with which it 
remained in equilibrium. 

For small departures from equilibrium 




a[c1 a 

< C >A 


A[C] 

TcT 


ML 

ML 


where [C] A is the carbon density in the atmosphere, [C]^ , that in the 

mixed layer, and £ ^ - 8*15; £ ^ , the so called "Revelle factor 11 , reaches 

the larger value at 0°C and the lower value at 30°C. The time needed to 

establish approximate CO^ equilibrium between the atmosphere and the 

mixed ocean layer is then a factor £ less than the 7 years needed to 

14 

attain equilibrium between their C contents; the latter essentially 
requires an exchange of the entire carbon reservoir of the mixed layer 
with that of the atmosphere. This time, 7£ yrs - 8 months is not short 
enough to maintain atmosphere-mixed layer CO^ equilibrium at all lati¬ 
tudes because horizontal surface motion exchanges water between hotter 
and colder regions in less than this time. However, it is short enough 
so that when averaged over the whole ocean mixed layer, CC^ equilibrium 
with the atmosphere is re-established after an atmospheric perturbation 
in less than a year. This is so very much less than the time for 
doubling the atmospheric CC^ burden from fossil fuel burning or the 
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residence time for a carbon atom in the mixed layer that we shall hence¬ 
forth assume atmosphere-mixed layer CO^ equilibrium is always achieved. 

The ocean can absorb more than the fraction £ - 10 ^ of any addi¬ 
tional CO^ injected into the atmosphere only because the rapidly respond 
ing mixed ocean layer is slowly replaced by fresher water which has not 
yet absorbed the added CO^ burden. Models for this replacement involve 
either or both of two mechanisms: 


i) The mixed layer exchanges water with layers directly 
below it, probably by the exchange of eddies — a 
process represented phenomenologically by a diffusion 
equation (box models). The time scale depends 
sensitively on the spatial separation between the 
mixed layer and the depth of the most relevant 
fresher layers. 

li) Cold water from the mixed layer can, in certain 

regions, drop relatively rapidly to much lower depths 

then spread horizonally and be replaced by a continued 

upwelling in the rest of the ocean. The largest such 

removal of surface water appears to occur in the 

Antarctic region of the Pacific — a downflow of about 
6 3 “l 

20 x 10 ms . This would be balanced by an average 

upwelling of A.A my/yr in the rest of the Pacific. In 

the Atlantic a similar upwelling might be expected 

6 3 “*1 

from a downflow of 10 x 10 ms in the Norwegian 
Sea. The time scale for replacing the mixed layer of 
thickness H * 75 m by this mechanism is 
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17 yrs 


( 2 ) 


75 m 


ML 4.4 m/yr 


We consider first only an estimate of the maximum consequences of 
mechanism ii) for ventilating the mixer layer. 

We define the total carbon contents of the atmosphere, mixed layer, 
and the deeper rest of the ocean by , and , respectively. 
Then, if CO^ is injected into the atmosphere-ocean at the rate F(t) 


f(t) 


C A + 


"ML 


+ C T 


(3) 


At present 



(4) 


with C of Equation (1). The factor £ will decrease if very much more 
CC >2 enters the system because of the near saturation of the mixed layer 
and the ratio C^/C^ decrease as increases. We shall approxi¬ 

mate the product 


» 


tC ML /C A 5 


(5) 
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by a constant equal to the present value « 0.7? * 10 ^ . The direct 

coupling between and Is model dependent because [C]^ can vary 

throughout the deeper layers. Therefore, the negative contribution to 
• 

Cp carried by the general upwelling accompanying the local downflow Into 
It from the mixed layer Is not simply related to . We maximize the 

effect of the upwelling from the cold water sinks by assuming 


°ML “ C ML ( " °°^ 


ML 


( 6 ) 


The reference time t » refers to the steady state which existed before 
fossil fuel burning significantly changed and thus , when C^O - ®) 

- 0 . This assumption that the upwelling water into the mixed layer Is so 
old that its carbon content is independent of recent increases from 

terrestrial surface activity during the past 150 yrs will, of course, 
exaggerate the effect of this mechanism in supplying surface water capable 
of maximally absorbing increasing atmospheric CO^. Because 

[V ] is very small, we can use Equations (4) and (5) 
in Equations (3) and (6) to obtain 


and 


C A - C A ( ~ 
T ML 


e' 


C A ■ C A ( " ~ ) 

f(t) =CL + 5')C. + — T - - C 

A ML 


(7) 


(8) 
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For the present exponentially increasing CC^ injection 


f “ f Q exp (t/x f ) 


and the solution of Equation (8) is 


(9) 


- C A ( - 


T f f o exp (t/T f> 


1 + C' 


I 1 + (1 + Aj 


-l 


( 10 ) 


With » 23 yrs corresponding to a 4.3% increase per year in CO^ injec¬ 
tion into the ocean-atmosphere reservoir, “ 17 yrs , and £' ■ 0.1 , 

Equation (10) gives 0.8 for the total fraction of injected CC^ which 
should still reside in the atmosphere. Because of the small £’, a t 

ML 

very much less than Is needed before would depart 

greatly from the value it would have when In equilibrium with a non¬ 
vent Hated mixed layer, represented by the first term on the RHS of 
Equation (10). This is particularly true in that the use of Equation (6) 
makes the whole discussion an overestimate of the cold water sinks. Thus, 
it is very unlikely that cold water sinks of presently estimated magnitude 
could by themselves account for the observations which show only about 
half of the CO 2 from fossil fuel burning still airborne. When, however, 
the continually accelerated rate of fossil fuel burning begins to 

diminish, T r will grow and the corrective factor Jl + £ f x(1 + £ f ) 

f Li ML 

can be important; but, if atmospheric CO^ has increased greatly by that 
time, the accompanying decrease in £' may no longer be small. 


39 


Digitized by LjOoq le 




B. A Pipe Model of the Oceans 


It is a sad commentary on the state of oceanography that we do not 
know of the applicability (if any) of two extreme models: 


(i) oceanwide diffusive and convective exchange processes 
adjoining vertical layers, and 

(ii) Processes at the ocean boundaries which alter the 
temperature and salinity (and hence density), with the 
newly formed water masses spreading into the ocean 
interior along isopycnal surfaces. 


The question enters the CO^ problem in a vital manner. The ocean sink 
associated with the formation of water penetrating to intermediate depths 
has a time constant of only 100 years; for bottom water this is 1000 
years. Since we are trying to make predictions for a century, this might 
be a vital difference. 

One-d models are rightfully in disrepute and we do not need to join 
the chorus of disclaimers (all of whom use such models). In fact, one of 
us (WM) is the originator of a particularly naive one-d model^, but in 
limiting this model to Antarctic bottom water and the resulting constant 
rate of upwelling by 1 cm day \ and ignoring the large intermediate sources, 
he now suspects that he has thrown out the baby with the bottom water. 
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1. A Pipe Model 


The ocean is modeled by a series of i pipes, from 1 * 1 at 
Che bottom to i - n representing the mixed layer at top. The tempera¬ 
tures, salinity, density and some other stuff (deuterium, CO^) in pipe 
i are designated T^ , , x^ . These are averages formed throughout 

the length of the pipe. Neighboring pipes are connected, with q^ de¬ 
signated the volume flux from pipe i into i + 1 . q Q * 0 (the bot¬ 
tom is taken as impermeable), and q Q " 0 (ignoring for the present 
evaporation and precipitation). 

Pipes are connected by conducting wires (Fig. 11-1), permitting 
the diffusion of heat (and salt and stuff), with the flux proportional to 
the temperature difference. This is to model turbulent diffusion in the 
ocean. Accordingly, the eddy coefficients are the same for heat, salt 
and stuff. 


3 

Finally, there is an import by m /sec of water with the 

i t i i • t 

properties T^ , , x^ . In general, T^ >* T^ and , 

t • 

Xi * x t ; but , for if it were otherwise, the inflowing water 

would be too heavy or too light and would sink to the appropriate depth 
in a few Vaisala periods. 

2. The Water Mass Makers 

In the present picture we take the view that the water is 
drawn from the few uppermost pipes (Q^ is negative for i near n) and 
then converted by "Water Mass Makers" and fed into the various lower 
pipes. What we are trying to model are some of the following processes: 
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a) The Norwegian Sea is filled with shallow Atlantic water, 
in winter the surface water is cooled by contact with the atmosphere, 
eventually the water column turns over forming a fairly homogeneous cold 
and moderately salty water mass which flows over a fill into a depth 
appropriate to its density; 

b) In the process of freezing at high latitudes, the newly 
formed sea ice is fresh, and the residual water is salty and cold. 

Some form of this process in the Weddell Sea is probably responsible for 
the formation of Antarctic Bottom Water; 

c) Extensive evaporation from the Mediterranean Sea results in 
increased salinity that leads to the formation of the densest waters in 
the world's ocean. In this discussion we use the oceanographic unit 

6 3 19 3 

1 Sverdrup * 10 m /sec * 3 x 10 cm /year 

3. Dilution 

3 

Suppose that on the average Q m /cm of Mediterranean water 
flows through Gibralter and down the continental slope into the Atlantic. 

In the process it mixes with the local water with the consequence that 

3 

by the time it arrives at equilibrium depth, the flow consists of Q v m /sec 
of diluted Mediterranean water, with Q* > Q . We will want, somehow, 
to include this dilution as part of the process of water mass formation. 
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4. 


Conservation 


We now get the following equation for various conservations: 


Mass: ^ + 


Heat 


: + vi + K (T i-i - V ' q i T i + K <T i ' T i+i ) 


Salt 


: q^jS^ + Q ± s i + K < s i_i " S i ) “ ^i S i + K (S i " S i+1 ) 


To this we add the equation of state: 


p i " f(S i’ V P i " f(S i’ V 


so that there are 4n equations in the 4n unknown q^ 1^, S^, 


In general, if we start by specifying and , and hence 

i 

, and then calculate the steady state , S , and hence, , 

t 

we will find that p^ f . We will need some procedures for satisfying 
the restraint p^ = p^ . 


5. Continuum 


,14 2 


We now divide the conservation equations by A * 3 x 10 m , the 
surface area of the world's oceans. (We should use A(z) . Use the no¬ 
tation q^/A = w^ (the velocity of upwelling). Set p * Q/A , the 
specific source function. Further, write 
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KAL 2 -1 

tc “ —r— m sec 
A 

so that ic has the dimensions of a diffusivity. We then write the dif¬ 
ferential equation 


-k + -t-(wt) - S 


dz 


2 dz' 


(ID 


for temperature, and similarly for salinity. Here 


AW 


i * Vi - M i' V A 


( 12 ) 


is the increment of upwelling velocity at depth z^ due to a source 
The source function 


n-1 

S " 2 AW i T l nPUt 6 ( z - z i) + W 1 T b 6(z) 
i 


W T_S(z-z ) 
n S s' 


and 


n-1 


W = 6(z)W- + AW. 0(z-z.) - W 0(z-z ) 
1 Smmd i i n 8 


(13) 


where 6(z) is the step function, and z q ■ 0 is the sea bottom, 

z “ z is the sea surface, 
n 8 
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The foregoing equation then becomes 


2 

d T dT 

-K + W ■ S - T dW 

,2 dz 

dz dz 


n-1 


y AW i < T i“ PUt “ T( Zl ))6(z- Zi ). 

(14) 


1-1 


We expect that the ocean can be modelled by a relatively small number 
of source depths, corresponding perhaps to the Antarctic Bottom Water, 
North Atlantic Deep Water and Mediterranean Water. Within any layer 
between two source depths we have simply 


d 2 T . „ dT 

-* 77 + w dl 

dz 


0 


(15) 


with simple analytical solution. 


The problem is to come up with reasonable source functions that 
lead to acceptable distribution T(z) , S(z) . The problem is mathematically 
trivial, but oceanographically very touchy. 
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c. 


Solutions to the Pipe Models 


Mathematically, the model described above may be set up as follows: 

we consider a single quantity, such as CO^ concentration, C(z,t) , which 

we take to be horizontally uniform and a function only of depth and time, 

in accord with our previous discussion. At the surface, z , we assume 

8 

C(z ,t) - C (t) , a given C0„ concentration in the mixed layer. At the 

8 8 Z 

bottom, z , we take C(z_,t) ■ C_ , a constant. In between the 
B Bo 

surface and the bottom are N sources of CO^, at depths z^ . . . z^ , 

injecting (X^ with concentrations C^(t) at rates AW^ . (AW^ has the 

dimensions of a velocity rather than volume per second because we normalize 

everything per unit horizontal area; AW^ is just the rate at which water 

is injected by the i C ^ "pipe". That is, AW^ = Q^/unit area. Thus 

between the i^ and (i+l) 8t "pipe" there is an upwelling of water with 

velocity , where AW^ ■ . To conserve water, these must 

then be an outflow of water at a rate AW^ * immediately below the 

surface z . This removes CO. with concentration C (t) . The geometry 
s L s' 

is illustrated in Fig. II-2.) C0 2 is thus distributed through the water 

column both by diffusion and by being pumped from the surface to various 
depths z^ at known rates, AW^ . The CO^ concentration then satisfies 
the diffusion equation 

2 n 

-K + ~ (WC) » - I 7 + 1 AW C 6 (z-z ) - W.C 6 (z-z ) (16) 

- Z oZ dt . .. I 1 1 1 8 8 

dZ 1*1 

with the Boundary condition 


C(z ,t) = C (t) (17a) 

S 8 

C(z B ,t) - C B (17b) 
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i 


C (t) c (t) 

s s 


z 

s 

W 1 

AW - - W. 
o 1 

z, 


„ _ 

1 

‘ w 

W 2 

AW X = W 1 * W 2 



C 

Z 2 

W 3 

aw 2 - w 2 - w 3 

t 

« 

W N 


Z N 


AW = W 

N N 


W NH " 0 


Figure II-2 

THE GEOMETRY ILLUSTRATING THE USE OF THE PIPE MODEL 
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For the case we are dealing with here, all the C ■ C (t) since the 

1 8 

pumping mechanism simply redistributes the surface CO^ concentration. He 
write (16) in a more general form since we will later apply it to tempera¬ 
ture and salinity, where, as described before, the pumping mechanism 
injects water with different temperature and salinity from the surface 
water. 


We note that C(z,t) is continuous across the layer boundaries, and 
the vertical derivative 3C/3z is not. In fact, from (16) we see that 


lim Kf-irrl 
e+0 


fei -fz y ( w i - w i + i)( c(z i* t> 

\ |z“Zj+e z«Zj-t J ' 


Thus we can also formulate the problem by writing 


C(z,t) - C^(t) 


z > z * z 
i-1 i 


where 


3 2 C 3C -3C 

1 . „ i 1 

_IC a 2 W i 3z 3t 

3z 


(18) 


with the Boundary conditions 


(i) 


C.(z ,t) - C (t) 

1 8 S 


(19a) 


(ii) 

(iii) 

(iv) 


°N+1 “ °B 


VV 0 - C 1+1 (Zj.t) 1 - 1 


< 1 %) 


(19c) 


■ < H > - W «)( C ‘U; e *) 1 - 1 -" 


(19d) 
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Equation (16) leads to a simple conservation law. By integrating 

from z„ to z we find 
B s 

A / c(z.t) dz - I - f I ) + (-V. + I 4 Vi) (20) 

z b \ i*. iv ' 1-1 ' 

If we have C. = C , i ■ 1 ... N , then the second term on the RHS of 
i s 

(20) disappears. This simply expresses the obvious fact that the rate of 
change of that CO^ in the water column equals the rate at which it 
diffuses in through the surface minus the rate at which it diffuses out 
into the bottom. 


Temperature, T(z,t) , and salinity, S(z,t) , satisfy the same 

equation with two modifications. First, we take the surface boundary 

conditions to be time independent so that T * T(z) and S * S(z) only; 

thus 3T/3t and 3S/3t disappear. Secondly, we assume that the pumping 

mechanism changes the surface temperature and salinity so that T i 4 T 

1 8 

and S g . Instead we invoke the equation of state and require the 

density p(T^,S^) of the injected water to equal the density 
^ptT(z^), S(z^)]^ of the ambient water. For given AW^, T^ and 
these conditions are used to determine the injection depths z^ 

The first step in solving an equation like (16) is to Fourier 
transform in time. We define 

00 

C(z,io) ■ J dte lu>t C(z,t) (21) 
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with inverse 


C(z,t) » j e -iu,t C(z,u>) 


( 22 ) 


Corresponding definitions obtain for C g (t) . Then (16) becomes, for the 


Fourier transforms. 


N 


-k (WC) - iti)C « -W.C 6(z-z ) + £ AW.C.6(z-z_) . (23) 

Z dZ Is S ,.11 1 


dz 


i-1 


Equivalently, (18) becomes 


» 2 c. 3C 

T2 + M liT -iuC i ' 0 

dZ 


(24) 


and the boundary conditions in (19) also hold for the Fourier transforms. 
The solution to (24) is 


C^(z,u>) 


A+( W ) 


a + (w)z 


(io)z 


A^(to) 


1 ... N+l 


(25) 


where 



(26) 


The boundary conditions in (19) then form a set of 2N + 2 linear 

equations which determine the A^(b>) in terms of C (u) and C_ . The 

i s B 

solution is then obtained by evaluating Eq. (22). For a large value of 
N , this is in practice a messy procedure, and we shall therefore, for 
purposes of orientation, first look at the case of a single source. 
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The model is shown below in Figure 11-3 


C (t) C 

s s 

Zg _ _ 

W 

W 

c 

s 


w 

W - 0 


z 


B 



Figure II-3 


For simplicity we chose z = and C_ m 0 . We also chose z.. ~ 0 . 

D D X 

In the two layers we then have 


Thus 


with 


ff I 

-kC 1 + WC X - iu>C^ “ 0 

ff 

-KC2 - iwC^ = 0 


W_ 

C x . ( A+ e“ + A_e'“) e 2k 
C 2 - (b + C 6z + B_e“ 62 ) 



B - f-\ 2 
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Since C 2 must vanish as z , we see B_ * 0 . (Note that 

/-Iio/k has a positive real part for u> in the upper half plane.) The 
remaining boundary conditions are 

az 


(v 


-az 


W 


e 8 + A e e^ K z 


8 C 8 (U>) 


A + + A_ 


ica(A + - A_) - (| + k/b) B + - WC g (u>) . 


(30a) 

(30b) 

(30c) 


These three equations can be solved for A + , A and B + 


Wz 

C. (z,<u) ■ e^ K C (to) 
1 8 


W 8inha(z 

i — i 

g -z) + ^ca coshaz + 

1 sinhaz 

Wz 

S 

” 2k 

l e 


+ <B^ sinhaz g + ko coshaz g 



Ola) 


C 8 (z,a,) “ T7w-r 


e^ Z C (to) 
s 


sinhaz + ica coshaz 
s s 


(31b) 


We also write the solution for the quantity tc 3C/3z| giitz , which is 


relevant for the conservation law in Eq. (20). (Note that k 3C/3z 


z=— 00 


vanishes.) 


3C 

3z 


Z“Z 


Wz 
_ 8 

¥ C (to) + Kae C (to) 
Is 8 


kci 8 inha z + 
( 8 


(i♦«») 


coshaz g ^e 


Wz 
_s 

2k 


- W 


[(y + <6^ sinhaz g + ko cosh<Z £ 


(32) 
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To make those somewhat formidable looking expressions more transparent, 
let us choose 


C (t) 
s 


C e 
s 


t/x 


0 


t < T 
t > T . 


Then 


C (to) 

s 


(itu + 1 /t)T 

r — _ 

s (io) + 1 /t) 


Since the solution C^(z,dj) and C^Cz.di) are proportional to C g (oi) , 
they both have a pole in the upper half oj plane at u> ■ i/t . This 
pole produces, on inverting the Fourier transforms through Eq. (22), the 
"steady state" contributions to C^(z,t) and C^(z,t) ; that is, the 
part proportional to . The coefficients of e t ^ T are simply 

residues of this pole. These residues are easy to compute, particularly 
in various limits. 


We can identify three times in the problem: (i) x , the growth rate 

of the mixed layer concentration; (ii) z g /W the time for transport of 

2 

water through the "pipe"; and, (iii) k/W , a diffusion time. 

In the limit where the diffusion time is very short — i.e., 

2 2 

t >> k/W and z /W >> k/W — it is easy to see that the "steady 
s 

state” solution is simply 
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The solution is simply constant in depth down to the source, and below 
that dies off exponentially. This is because in the upper layer the 
CC >2 is uniformly mixed by the "pipe", while in the lower it spreads only 
by diffusion. 

The opposite limit, in which W -»■ 0 , of course yields simply 

z-z 
_ s 

(z,t) - C 2 (z,t) - ; (34) 

now, the only spreading is by diffusion everywhere. 

These steady state solutions, however, are relevant physically only 
if the time for water transfer through the "pipe" is short compared to 
the growth time t of the mixed layer concentration; i.e., if z g /W << t . 
For the real ocean, however, we have z g ~ kilometers, W ~ a few meters per 
year, so z /W ~ hundreds of years, while the characteristic time for CO.. 

8 Z 

growth is only about 20 years. Thus the pumping mechanism cannot in 
fact keep up with the CO^ growth, and hence the steady state situation 
is Irrelevant; we must look for transients. This means we cannot confine 
our attention only to the pole as W ■ i/x ; we must look more widely at 
the u dependence. 
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This is not difficult to do if we are willing to assume that the 
"pipe" transport dominates diffusion — that is, if k is very small. 

It is easier to look directly at Eq. (32). When k -*■ 0 , we find from 
(32) that 


K 


3C 

dz 


z=z 

s 



(35) 


+ terms which vanish with k 


Upon inverting the Fourier transform, and invoking the conservation law 
(20), we find 



(36) 


Thus, if S/W << t , the right hand side is just 


s Tt c s (t) 


in conformity with the steady state solution in which C(z,t) is 
constant in z in the upper layer of the ocean. But if S/W >> x , 
then the right hand side is 


WC (t) . 

s 


Thus, 


in the physically relevant case we find, if diffusion is small 

z 

s 



C (z,t) dz - WC (t) 
s s 


( 37 ) 
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This case has been analyzed by Ruderman (see § II-A of this report — 

Eq. (37) above is Ruderman"s Eq. (6)), so we do not need to study it 

further here. Numerically, for this case, the effect of the "pipe" is 

not large. For it to be large, one would have to be in the regime where 

z /W << x instead. 

8 


Discussion 


A full exploration of the "pipe mode" for ocean circulation will 

require detailed oceanographic data on the movement of water masses, 

particularly those generated in the cold regimes. Estimates of a down- 
6 3 —I 

flow of 20 X 10 m S in the Antarctic region of the Pacific and 
6 3 —1 

10 x 10 m S from the Norwegian Sea lead to the conclusion that the 
cold water sinks cannot account for the observation that only half of 
the carbon introduced into the atmosphere by burning carbon based fuels 
remains in the atmosphere. If mass flow from the surface downward to 
substantial depths is to provide a significant sink, then either the 
Antarctic and Norwegian flows are underestimated or other areas of the 
oceans not yet identified contributed to the removal of carbon from the 
atmosphere. 

Within the uncertainties of the "pipe model", type calculations, it 
appears that the conclusion that the ocean is a limited sink for carbon, 
as derived from one dimensional turbulent diffusion calculation is 
supported. It is difficult to understand how the oceans take up one half 
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of the fuel generated carbon. A requirement that oceans provide a sink 
for additional carbon derived from the biosphere would raise further 
difficulties in understanding both models of ocean processes and the 
oceanographic data that have been used in constructing these models. 
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Ill 


RADIATIVE LIMITS ON CLIMATE 


The sun's radiation in the visible part o£ the spectrum, except for 
that fraction which is directly reflected back out into space, warms the 
earth. For an earth in radiative equilibrium, the solar radiation which 
interacts with the atmosphere and the earth's surface must be balanced by 
an equal amount of outgoing infrared radiation. The temperature of the 
earth as observed from space, the "effective temperature", is determined 
by the amount of energy which the earth must lose to remain in radiative 
equilibrium. The effective temperature is lower than the surface tempera¬ 
ture since all parts of the atmosphere radiate into space as well as the 
land and oceans. The actual surface temperature and the temperature pro¬ 
file of the atmosphere are determined by the interplay of radiative pro¬ 
cesses, convective transfer of heat and the heat released in the phase 
transitions of water. Altering the radiative properties through changing 
the concentrations of infrared absorbers will change climate. 

In this section, we explore two models of the atmosphere is which 
only the radiative processes are taken into account. These models are a 
first step in a hierarchy of atmospheric models of increasing complexity. 

The study of the effect of carbon dioxide of the radiative properties 
of the atmosphere has a long history tracing back to the work of Tyndale 
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in 1861. More recently a number of studies have explored the impact of 
changing carbon dioxide content on climate but the emphasis has been on 
the 15 ym absorption band of (X^ (see for example Plass^). In this re¬ 
port we emphasize that the 16 ym region is almost opaque at current con¬ 
centrations of C (>2 an d that the bands in the 10 ym region are of greatest 
interest. Furthermore, other minor constituents of the atmosphere enhance 
the warming effect of CC^ by blocking outgoing infrared radiations. 


A. Gray Homogeneous Model of the Atmosphere 

The troposphere is not in radiative equilibrium nor is it even appro¬ 
ximately gray (meaning that the absorption coefficient is Independent of 
frequency). Further, it is not vertically homogeneous, since I^O is dis¬ 
tributed with a scale height of about 2 km, C (>2 with 8 km, and 0^ existing 
mainly in the stratosphere. Hence it is surprising that a gray, radiating 
homogeneous model has any value at all. Certainly this kind of model would 
not yield a very useful temperature-height profile. 

Such a simplified model can be useful, however, for exploring the 
importance of changes in the greenhouse effect caused by changes in the 
composition in global heat balance models of climate. 
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1. The Gray Model 


In the gray model only two temperatures are important: The 
effective emission temperature, T^ , of the planet is fixed by a 
balance between the Incident solar flux and the planetary thermal emis¬ 
sion. For a rotating planet of radius R. 

4irR 2 crT e 4 - (l-A)nR 2 (irF 0 ) (1) 


where A is the omnidirectional albedo and ttFq is the solar flux out¬ 
side the atmosphere. The surface temperature, T , which is assumed to 

8 

be the immediate source of radiation as far as atmospheric heating is con- 

2 

cemed, is related to T @ in a gray model by (Chamberlain, 1978, p. 11) 

T s ■ T . « + ! V • (2) 


where t is the vertical opacity of the atmosphere. 

s 


If the ground emits isotropically with intensity B. (T ) , the 

A 8 

flux transmission directly to cold space is 


2 IT 1 


oT’ Q(t) 


■fdX B j( I s )fdij 


Mdy 


o o 


- oT 4 2E (t) « oT 4 e 
8 3 s 


(3) 
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where E^ CO is the exponential integral, which we approximate with a 
simple exponential. As x varies from <*> to 0 and the transmission 
function Q * 2E2 (t) varies from 0 to 1, the coefficient g varies from 

1.2 to 2.0. For HjO and 0^ we choose 6 ■ 1.66 , which gives the 
exact value of E^ when Q - 1/2; for C0 2 at 15 pm , we take 0 - 1.2 

For a realistic atmosphere the same kind of considerations 
hold except that T is a function of X , and 


n 


Q 




(4) 


where the summation is over narrow wavelength intervals and p^ is the 
fractional flux in the i interval. 




oT 


(5) 


2. Earth’s Infrared Absorption 

Between 0.54 and 8 pm 1^0 has strong vibration-rotation bands 
and longward of 15 pm occur rotational lines from the ground vibrational 
state. The Earth's thermal radiation is mainly confined to the 5 to 


64 


Digitized by 



30 pm region. The H 2 0 molecule is thus a strong absorber over the 
thermal spectrum, except for the important 8 to 18 interval, 

where the Earth’s emmision peaks (at 290°K, X *10 pm). The 12 
to 18 ym region is dominated by the CC >2 v 2 “ pending) mode fundamental, 
leaving an 8 to 12 ym transparency window. A portion of this window 
is blocked by the 9.6 ym band of -»■ 0^. Greenhouse heating of the 
atmosphere is thus dependent on composition in several distinct ways. 

First, water vapor is dominant over a wide wavelength interval, 
although in much of this region it is not optically thick. Hence 
atmospheric cooling to space is critically dependent on ^0 abundance 
and even the vapor in the stratosphere is significant. (The latitudinal 
variation of ^0 is the main cause of the difference in temperature 
profiles between polar and tropical zones). 

Second, COj fills a large part of the ^0 window, and current 
C (>2 abundances are enough to make the 15 ym region optically thick. 

Major CC >2 abundance changes will alter the greenhouse heating mainly 
by C (>2 absorption in the 10 ym bands. 

Finally, minor constituents that absorb strongly in the 8 - 
12 ym region can be critically important. Nitrous oxide (^0) and a 
wide array of hydrocarbons are potentially important, and the gray 
model can be used to estimate their effects. 
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The abundance of water vapor typically varies from 0.2 to 

2 

2gm/cm (or cm of precipitable water). We adopt a global average of 
2 

lgm/cm . With sufficient accuracy H^O may be regarded as completely 
opaque except in the near infrared and visual (P^ - .03) ; and 8.5 to 
12.5 where T « 0.10 and p = 0.29 ; and 12.5 - 15 pm where T < 1 and 
p = 0.14 . However, in the latter region H 2 0 does not contribute to 
atmospheric heating since it is dominated by 00^, and the H^0 therefore 
can be neglected from 12.5 to 18 pm with p * 0.26 (Kondratyev, 1969, 
p. 119). Thus Eq. (4) gives the absorption function due to water alone 
as 


A(H 2 0) * 1 - Q(H 2 0) = 0.45 (6) 

As C0 2 increases the mean opacity within a fixed wavelength 
interval increases, but also the band effectively widens as weaker ab¬ 
sorptions become stronger. Here we will not consider the broadening of 
the 15 pm band because a more important optically thin effect is the 
emission by the 10 pm bands. Thus the treatment of 15 pm + 10 pm together 
should give a good representation of (X> 2 . 

For C0 2 a volume mixing ratio of f = 333 ppm corresponds to 
an equivalent thickness at STP conditions of 

£(C0 2 ) = 266 atm-cm (7) 
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An empirical formula for the mean transmission in the nearly 
saturated 15 pm band is (Kondratyev, 1969, p 126)^ 



exp(-B 'a^ b ) 


( 8 ) 


where the absorption data near £ “ 300 atm—cm yield, where $’ ■ 1.2 , 
a = 0.361 , and b ■ 0.23 . Then the total absorption due to <X> 2 at 
15 ym (or more precisely between 12.5 and 18.2 ym) is 


A(C0 2 ) - 1-Q(C0 2 ) > p(15y) [1-exp(-3' a£ b )] - 0.27(0.791) = 0.214 

(9) 


For the weak C0 2 band absorption we use the treatment in 
Section V, below, for minor constituents. The properties of the two 
intercombination bands at 10 ym are listed in Table X1I-1 and from 
Eq. (32), we have A(10y) - 0.025. 

For 5 * 0*3 atm-cm of ozone in the stratosphere we find 
3 

(Kondratyev, 1969, p 137) over the band (9.4 to 9.9 y) a mean transmission 
of 


or 


^Q(9.6y)^ = exp(-2.8S) 


( 10 ) 


A(0 3 ) - 1 - Q(0 3 ) - p(9.6y)(1 - e -2 * 8 ^) - 0.035(.57) - 0.020 (11) 
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Table III-l 


ABSORBERS IN 8 - 12 pm WINDOW 


Species 

Band, A(pm) 

^0 Nitrous oxide 

v x 7.8 

CH. Methane 

4 

v. 7.66 

4 

NH^ Ammonia 

v 2 10.5 

CF 2 C£ 2 F-12 

v x 9.1 

v, 8.7 

D 

v 8 10 ‘ 9 

CFCi 3 F-ll 

v x 9.2 

v 11.8 

4 

Hydrocarbons 

— 

co 2 

2v 2 ~v 3 9.4 

10.4 


f (PPm) 

S(cm) 

NS(cm" 1 

0.28 

O 

1 

•^1 

60. 

1.6 

7xl0 -18 

240. 

6xl0 -3 

-10‘ 2 °(?) 

-1. 

-4 

1x10 

5X10' 17 

0.1 


3xl0 -17 

0.06 


5xl0 _17 

0.1 

•4 

1x10 

3xl0~ 17 

0.06 


6xl0~ 17 

0.1 

< 4xl0 _3 

<5xl0 _17 

<4- 


-22 

8.6x10 

6. 

333 

-21 

2.7x10 

19. 


Note: For uniform mixing the Integrated abundance is N » fN ^ = 

atm 

f x 2.15 x 10 25 cm" 2 . Since f = N/N - N/2.69 x 10 19 and S.T.P. 

o 

thickness is £ = N/N * fN _ /N , we have 

o atm o 

f(ppm) = 1.25 5(atm-cm) . 

If number in last column exceeds 10, a doubling of present abundance 
would produce AT > 1°C . 
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If the various substances do not overlap seriously In their 


spectral regions of absorption, the A's can be added. Thus from Eqs. (4), 
(6), (9), and (11), we have (with 3 ** 1.66) 

-8t 

Q net - e 8 - 1 - [A(H 2 0) + A(C0 2 , 15m) + A(C0 2 , 10y) + A((> 3 )1 


giving 


t - 0.748 . (13) 

g 


- 0.29 
( 12 ) 


3. Earth*8 Emission and Surface Temperatures 

With an albedo A • 29% , Eq. (1) gives T g * 257°K . Then 
Eqs. (2) and (13) yield 

T - 287.2 (14) 

8 

which is close to the world average (~290°K). It is probably fortuitous 
that the agreement is so good, although it does suggest that the 
gray model may provide some insight into climate changes due to changes 
in t 

g 
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4. Variation of Ground Temperature with CO? Abundance 


Differentiating Eq. (12) gives, with Eqs. (9) and (32), below. 


8e 


-g T 

dx 


p(15n) exp(-3'a5 b ) B'abC b 


+ 1.08 x io “ 3 S-N . 

v o dx 


(15) 


In terms of the mixing ratio or volume concentration 


f — 
df 


[co i dx , p(15y)g , ab exp(-8'aS b K b 08x10 S-N Q g 
[ 2 ' 0e _6T 8e' 6T 


4.21 x IQ -2 + 7.5 x 10 -2 - 0.117 


(16) 


Similarly from Eq. (2) we obtain 


dT 


3T 


3T 


dT 16(1 + .75 x) 3/4 16(1 + .75x) 


34.5 


(17) 


and therefore 


dT 
_s 

df 


4.0°K 


(18) 
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If the CO content were doubled, the gray model would give T = 290.0 

S 

or AT g = 2.8°K . This value compares nicely to values obtained from more 

elaborate models. In reviewing a number of recent studies, Schneider 
4 

(1975) concluded that a doubling of CO, would produce a AT of 1.5 to 

fc s 

3°K . As C0^ is increased further, the temperature will continue to rise; 
although the 15 pm bands become more saturated, the 10 pm bands are 
effectively closing the 8 - 12 pm window, as is illustrated in Fig. III-A. 

The energy-budget climate model (§ IVA) relates the infrared 

radiation to surface temperature, T (x) , by an empirical relation 

s 

I i oT g (x,f) - A + B[T g (x,f) - 273] (19) 

where x is the sine of the latitude. If the solar constant is indeed 
constant but f is allowed to vary, then < T g > (averaged over the 
planet) remains constant with f . Representing the latitudinal varia¬ 
tion of surface temperature by 

T (x,f) - T (f) + T,(f)P,(x) (20) 

8 O Z a 


we have 


<T g (f^> - T Q (f) (21) 

Even though A and B are each independent of local surface temperature 
in accord with Eq. (29), they must nevertheless vary with the global mean 
^T g (x,f)^ in order to preserve a constant ^T g ^ . Thus 
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♦ 


* 




*• 


4 


Figure III—A 

TEMPERATURE VERSUS C0 2 ABUNDANCE IN A GRAY MODEL 
IS IN GLOBAL RADIATIVE EQUILIBRIUM 


The inflection point around 750 - 1000 ppm develops because the 15 Mm band is becoming 
saturated but the 10 Mm bands are not yet strong enough to produce appreciable heating. 



« 


4 



( 22 ) 


A(f) + B(f) [T (f) - 273] - const, 
o 


As a practical method of computing a global model with a modi¬ 
fied f , we note from Eq. (2) that a gray radiative model valid at every 
latitude would yield a latitudinal dependence of flux of 


I[f,T g (x)] - I|f,T g o) j 


oT 


<o)4 


1+.75 T 


g 



+ . . • 


(23) 


where T ^ is a convenient reference temperature. For T = 273 , Eq. (23) 
8 8 

can be compared with Eq. (19) to obtain the A and B dependence on T 

g 

or f . This result suggests that A and B are scaled by 


A(f) 


A(f )(1+.75 t ) 

° _8i° 

1+.75 t 

g 


(24) 


and 


B(f) 


B(f )(1+.75 T ) 
o _ g»° 

1+.75 t 

g 


where T corresponds to abundance f 
g.o o 


(25) 
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when Che infra- 


The mean surface temperature of the planet T q , 

red flux is a specified function of latitude, Eq. (19), will not be the 

same as the simpler radiative mean, Eq. (2). If the optical thickness is 

modified by Ax , the radiative change is given by Eq. (17), whereas the 
8 

latitude dependence gives, by Eqs. (22), (24), and (25), 


AT 
_c 

At 


g 


3 /A/B + T o -273\ 

4 ( 1+.75 , ) 


(26) 


2 2 

With A ■ 211.2 W/m and B ■ 1.55 W/m °C , Eq. (25) is about twice the 
rate of Eq. (17). 


5. Radiative Effects of Minor Constituents 

The addition to the atmosphere of a minor constituent that ab¬ 
sorbs in the 8 to 12 pm window could be important. In this region 

4 

ttB^/oT^ ** .060 - .070/p and hence Eq. (5) is 

P ± « 0.065 AA y , (27) 

where AA^ is the width of the absorption in microns. We will estimate the 
minimum amount of a trace gas that could produce a 1°K change in global 
temperature. 


For an optically thin amount of absorber, we have simply 


A 


i 




^Pi< T i> » 


(28) 
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where <r> is the mean value in the band, or 


which gives 



dA 


(29) 


A 


8(.065) SI 



dA 


- .108 N S 

A 


(30) 


where N is the integrated column density (molecule/cm ) and the 

2 

integrated absorption coefficient or strength (in units of cm micron). 
The more common units for strengths are cm (cross section in cm; 
spectrum in cm "S: 


S~ 


v 



dv 


(31) 


Hence at 10 pm , 


A 


0.108 SI S~ 

v 



- 1.08 x 10 -3 SI S~ 

V 


(32) 


The maximum band strengths are not likely to exceed 5 x 10 ^ cm (see 
Appendix A, Ruderman). 
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From Eq. (12) we have 



-Bt 

Be ® 


1.66 (0.15) 
0.25 


(33) 


and with Eq. (17) the T g - A relation becomes 


dT 
_s 

dA 



dA/dx 

g 


= 138 


(34) 


Thus for 1°K changes in T g , we are searching for increments in A the 

_3 

order of 7 x 10 , or from Eq. (32), 

AW * - 7 X1 ° - 3 - 17 - 1*3 x 10 17 cm" 2 (35) 

1.08 x lO"'* x 5 x 10 

25 -2 

The integrated column density of the atmosphere is 2.15 x 10 cm 

so we are concerned with mixing ratios of minor constituents greater than 


Af ~ 5 x lO -9 (36) 

or a few parts per billion (ppb) for substances that have a single 
strong absorbing band in the 8 to 12 pm window. 


The halogenated methanes (Freons) have been cited as a striking 
example of the effect we are examining here (Ramanathan 1975)."* The 
two Freons F-ll and F-12 together have five bands with strengths 
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S -» 3 x 10 ^ cm . The 1975 abundance was 0.1 ppb; an Increase of 10 
s 

times would place these substances In the climate critical category. 

Most substances will not have bands In the window with strengths 
close to the maximum. From Eq. (35) the pertinent quantity for bands in 
the window is 


AW • S > 10 (37) 

if an absorber is to produce a 1°C change on future climate (see Table 111-1). 
Both N^O and CH^ sit in regions where there is already strong absorption 
and they have to be analyzed (as we did CC^) allowing for the present 
near-saturation. 

All of the substances in the Table except the hydrocarbons 
have been examined with a more elaborate radiative model (Wang, et al, 
1976)^and found to be marginally important, if their abundances increase. 

For N 2 O, CH^, and NH^, factors of two are important. 

More attention needs to be directed toward the global abundance 
and long term growth rate of hydrocarbons, which are among the primary 
urban pollutants, and whose emission is largely uncontrolled. The abun¬ 
dance in Table III-l applies to the total hydrocarbon population; the strength 
listed is the maximum from Appendix A and probably is much smaller. 

Aldehydes (a hydrocarbon family) are partly destroyed near their urban 
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source, where they participate in the production of oxidants. Away from 
high concentrations of NO^, the aldehydes disappear by photodissociation 
and attachment to aerosols. 

6. Destruction of Trace I.R. Absorbers 

Many of the hydrocarbons are not only being emitted at increasing 
rates, but their natural sinks are becoming saturated. The principal 
sink for methane is hydroxyl: 

CH^ + OH -*• CH 3 + H 2 0 . (38) 

Indeed, OH plays a crucial role in scavenging a number of the less 
soluable trace gases, such as CO, H 2 S, S0 2> and the partially halogenated 
methanes, CH x and CH^r^. A depletion of OH may thus aggravate a 

build up of infrared absorbing pollutants and their inpact on the climate. 
At present there is a danger of such a depletion occurring because of the 
increasing production of CO, which is removed mainly by 

OH + CO -*• C0 2 + H (39) 

9 

The present CO production from fossil fuel burning is 0.4 x 10 ton/year 

9 

whereas the natural production is 1 to 3 x 10 ton/year. The problem 
is amplified by the fact that the main natural source of CO is the oxida¬ 
tion of CH^. Following reaction (38), oxidation proceeds by 
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(40) 


CH 3 + 0 2 + H + CH 3 0 2 + M 


ch 3 o 2 + no -*■ ch 3 o + no 2 


(41) 


CH 3 0 + 0 2 -► H0 2 + CH 2 0 


(42) 


Finally, formaldehyde is dissociated into H 2 + CO, or it may condense 
onto aerosols and ultimately rain out. 

The principal source of tropospheric OH is photodissociation 
of tropospheric ozone: 


0 3 + h v (A < 4110A) -*• 0( D) + 0 2 


(43) 


followed by 


H 2 0 + o( D) -► 2 OH 


(44) 


# 






Thus present levels of man-made CO emission are large enough 
to perturb the C0-0H-CH^ cycle. Reactions that recycle OH, such as H0 2 + 
NO — OH + N0 2 » may account for some 30 per cent of the present OH. 

7. Feedbacks: Comments on Water and the Stratosphere 

Earth is the water planet and water affects climate differently 
in its different manifestations: vapor, clouds, sea, sea-ice, land-ice and 
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snow, and foliage. Common to all models is the acknowledgement of the 
effect and inability to do anything useful about it. It is not unlikely 
that the response of water to a change in CO 2 is Che dominant factor in 
climate change. If so, we do not even know the algebraic sign of the 
change. 


Also physical processes in the stratosphere can feed back on 
surface temperatures. The water abundance in the stratosphere is only 
several parts per million and is probably fixed by the vapor pressure at 
the tropical tropopause [although HgO can be made in the stratosphere by 
reaction (38)]. A change in the content of water due to a change in 
stratospheric temperature will affect the earth's greenhouse. An altera¬ 
tion of the ozone abundance will change the absorption at 9.6 lim, which 
is the main absorber in the 8 - 12 pm window. Thus an increase in tro¬ 
pospheric CH^ could affect stratospheric HO^ chemistry, the ozone abun¬ 
dance, and thence the climate. 


B. A Simple Band Model for Infrared Emission from the Terrestrial 
Atmosphere 

The energy budget climate model described elsewhere in this report 
requires knowledge of the infrared flux (I) radiated to space by which 
the Earth balances incoming solar radiation. The flux I must be given 
as a function of the local surface temperature (T g ) and a few other at¬ 
mospheric parameters such as the concentration, lapse rate, etc. 
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A proper treatment of Infrared emission from the Earth's atmosphere in¬ 
volves a multitude of complex factors, viz. Paltridge and Plott (1976) 
or Kondratyev (1969). Our hope here is to develop a model simple enough 
to be tractable yet incorporating enough observational and theoretical 
knowledge of the atmosphere so as to be responsive to changes in the at¬ 
mosphere such as a changing CO^ content. While this effort is intended 
primarily as an input to the aforementioned climate model, several in¬ 
teresting features are revealed by exercising the model. 

Bather than letting radiative transfer determine the temperature 
structure of the atmosphere we will assume that the atmospheric structure 
is known and simply let this assumed structure radiate to space. Our 
assumed structure is based on the U.S. Standard Atmosphere, 1976. In 
this sense the model can be seen as a perturbation on this Standard 
Atmosphere. 

Although we compare model results, i.e. values of the infrared (IR) 
flux radiated to space, with satellite observations and other empirical 
data, the model has not been "turned" to observational data except by the 
introduction of a normalizing constant for the total infrared flux and 
the use of empirical values of atmospheric parameters, e.g., lapse rates, 
scale heights, surface densities, etc. Where possible, these parameters 
are taken from the U.S. Standard Atmosphere (1976). 
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The model is briefly summarized In terms of three assumptions: 

• The atmosphere Is locally plane parallel with a tempera- . 
ture profile T(z) given by a piecewise linear curve con¬ 
sistent with the U.S. Standard Atmosphere (1976). For 
surface temperatures T g other than 288°K, the lapse rate 
is held constant and the tropopause temperature held con¬ 
stant at 217°K (see Fig. III-l). Structure above the tropo¬ 
pause is also held constant. 

• The atmosphere radiates 1R emission only in the 8-100y 
(100 - 1250 cm 1). This spectral range is divided into 
nine bands chosen so as to reflect salient features of 
the terrestrial IR spectrum and their changes over dif¬ 
ferent parts of the Earth. The bands correspond to the 
atmospheric window and the absorption bands of H^O, CO^ 
and 0 3 (see Table III-2 and Fig. III-2). 

• The IR emission of the Earth's surface arises from a 

black body at the surface temperature (T ) . Emission 

s 

in the optically thick CC >2 and 1^0 absorption bands 
arises from an atmospheric layer whose temperature is 
appropriate to an optical depth of 2/3 measured ver¬ 
tically downward from the top of the atmosphere 
(Chamberlain, 1978). The 0^ band and optically thin 
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Figure 111-1 

SCHEMATIC DIAGRAM OF TEMPERATURE PROFILE USED IN NINE-BAND IR 
EMISSION MODEL FOR THE TERRESTRIAL ATMOSPHERE 

The ordinates are height (z) or optical depth (r). Optical depth r is measured down¬ 
ward from the top of the atmosphere. The plot of T(z) is taken from the U.S. Standard 
Atmosphere, 1976. The dotted lines indicate perturbations for surface temperatures 
above and below 288° K. The r =2/3 levels for various bands are indicated for the case 
of T s =288° K, relative humidity = 50%, CO 2 content 9 332 ppm. 
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BAND CHARACTERISTICS 
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log AX 



- 1/X (in cm" 1 ) 


• Figure 111-2 

INFRARED BAND ABSORPTION OF ATMOSPHERIC GASES 


• The IR emission band (7T= 200-1000cm” 1 ) in which the model is allowed to radiate 

is shown as being composed of nine bands: atmospheric window (6), Og (8), COg 
(3,4, 5, 7, 8) and HgO (1,2,9). A black body curve for 250° K is shown normalized 
to the upper boundary of the graph. The plot is modified from that given by Allen 
(1976). 
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C0„ bands are handled as thin absorbing layers over the 
surface which radiates as a black body at T g . 

1. Model Description 

a. Division Into Bands 

We divide the terrestrial IR radiation spectrum into nine 
bands as shown in Figure III-2. This appears to be about the minimum number 
one can use and still hope to include the major effects of the absorbing 
gases. In Table III-2 we have summarized the characteristics of these 
bands. In band 1 we include the strong (50p) absorption band of water 
vapor. Band 2 also contains water vapor absorption, but considerably 
less strong than band 1. The strong CC> 2 , v 2 bending mode band at 15y 
is accounted for in bands 3, 4 and 5. 

The 8-12p atmospheric window is divided into four bands. 

Band 6 is completely transparent and thus radiates at an effective tem¬ 
perature equal to the surface temperature T . One can consider radiation 

s 

in this band to be from a grey body with emissivity e . Since e ~ 0.95 , we 
shall simply take it as unity in our calculations here. Band 7 presently 
contains only weak water vapor and CO^ absorption. We have added this 
band in anticipation of a major rise in atmospheric C0 2 content. As 
noted by Chamberlain elsewhere in this report, the major effect of increases 
in C0 2 content, so far as the greenhouse effect is concerned, will proba¬ 
bly lie in making currently transparent portions of the atmospheric window 
more opaque, rather than making the 15y C0^ band even more opaque. Band 8 
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contains the 9.6y ozone and 9.4y COabsorption bands. Band 9 contains 
relatively weak water vapor absorption. 

b. Total IR Radiation to Space 

The total IR radiation to space is given by a summation 

over 9 terms: 


9 

F - a£ IIB-Iv^ T^,...)^ (45) 

1=1 

where the Planck function for V in cm~\ T in K is 

B~(V,T) * Cl V 3 /[exp(c 2 \)/T) - 1] Wm ^sr 3 (cm 3 ) 

c - 1.911 x 10" 8 
c 2 = 1.439 

and 

V ± ■ band center wavenumber cm 3 
Av^ = bandwidth in cm * 


The factor ir arises when converting the radiance or specific intensity 
B~ to a flux through a plane parallel atmosphere. The constant a , of 
order unity, is used, as discussed below, to normalize the value of I to 
some desired value for a specified set of conditions. For example, it 
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could account for IR outside the 

temperatures are determined 

T.(T ,...) Is a function of the 
1 s 

atmospheric parameters, e.q. CO 


100-1250 cm ^ band. The effective 

for each band as discussed below. 

surface temperature T as well as 

s 

content, relative humidity, etc. 


Each 

other 


2. Effective Temperature Calculations 

The effective temperatures for optically thick bands are 
determined by noting the altitude (or equivalent pressure level) at 
which the mean optical depth measured vertical downward from the top of 
the atmosphere in band i is 2/3. The temperature T^ at this level is 
then the effective temperature at which the i C ^ band radiates. The 
value 2/3 for the optical depth of the effective radiating level is 
derived by Chamberlain (1978). 

a. Transparent Atmospheric Window (800-925 cm \ Band 6) 

In this completely transparent portion of the atmospheric 
window the effective radiating level is at the surface (z^ * 0) and the 
effective temperature simply the surface temperature T^ * T g . 


b. Water Vapor Bands (1, 2 and 9) 

For altitudes below about 30 km the lines in the water 
vapor bands of interest here are not well separated and overlap consider¬ 
ably. Goody (1964) has developed a random model to account for this over¬ 
lap in which the mean transmission T in a spectral range Av is given 
by 
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where is the total equivalent width of the lines in the range Av . A 

mean vertical optical depth of 2/3 implies a mean transmission T * 1 
-exp(-2/3) * 0.487 above the effective radiating level, z^ . Thus z^ is 
defined by making such that 

V^/AV - - In (T - 0.487) = 0.72 (47) 

is calculated differently depending on whether or not the lines in 
the Interval Av^ are weak (unsaturated) or strong (saturated such that 
additional absorption can only occur in the wings of a line) - see 
Houghton (1977) for reference. 


Following Houghton (1977, Ch. 4) the equivalent width for 
th 


weak lines over the i band is given by 


W i = ^ W j J S j C P dzRsS i / c.pdz 

^ ^path path 


(48) 


where is the equivalent width over the j*"* 1 subinterval inside the inter¬ 

val Av , it is the path length, c the mass mixing ratio of the absorber 
and p the atmospheric mass density. Here we have assumed s^ constant 
along the path. The quantity s, is the line strength for lines in the 
j t * 1 subinterval. ^ s^ in units of cm *(kg m ^) 

j 
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For strong lines, the equivalent width is 


j j 


J* cpdz 

~’E ( v / /2 

J' cpdz 

path 

j 

path 


1/2 


where we assume 8 jYj is constant along the path. 


Taking account of pressure broadening by letting the 


half width Yj * PY Q j/P 0 we have 


■ *rf’/ 

j [path 


(8 i Y oj )(cPp/p 0 ] ! dz 


r-•■■[/. 


(cPp/p )dz 
o 


2 


1/2 


(49) 


where p is the local pressure, p Q the standard pressure (1.013 x 10 Pa), 

E l/2 -1 -2 -1 

(SjY„ 4 ) in units of cm (kg m ) , and y„ 4 is the half width 


'oj‘ 


oj 


j 


, th 


over the j subinterval at STP. The quantities S and R are conveniently 
tabulated (in cgs, not mks unitsl) by Houghton (1977, Appendix 10) for water 
vapor, carbon dioxide and ozone. The values of S or R , appropriate 
to the effective temperature and summed over bands of interest, are shown 
in Table 1. S and R are temperature dependent and especially for weak 
lines this dependence can be important. To include the temperature varia¬ 
tion we have used linear Interpolation between the S and R values given 
for 220, 260 and 300 K by Houghton. This linear interpolation when com¬ 
bined with a simple iterative computational scheme yields values of R and 
S which are consistent with the effective temperatures determined. 
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It now remains to calculate the integral quantities in 

Eqs. (48) and (49) for both weak and strong lines. In the case of water vapour 

we let cp = p = p exp (-z/H ) with a scale height H of about 2 km. The 
w wo w w 

surface mass density of water vapour is discussed below. For weak lines 

from Eq. (48) we have 


W 


i 


S ,P 
i wo 



path 



S,P H e 
i wo w 


-z/H 
i w 


(50) 


where the path is from the effective radiating level z^ vertically up¬ 
ward to infinity. 


For strong lines from Eq. (49), letting cp - p , as above, 

w 

and (p/p Q ) ■ exp(-z/H) where H is the scale height for the atmosphere as 
a whole, we have: 



(51) 


3 3 

Using H * 2x10 and H = 8.5x10 m we have H H/(H +H) = 1619 and (51) 
w w w 

becomes 


1/2 


W.-80.5 R p 
i i wo 


exp(-z i /3238) 


(52) 


91 


Digitized by v^ooQle 



We must now estimate the surface mass density of water 

vapour p . A number of workers (e.g. Manabe and Wetherald, 1975) have 
wo 

used a constant relative humidity formulation. Following this sugges¬ 
tion we shall calculate p so as to maintain a constant relative 

wo 

humidity r (r ** 0.5 for 50% relative humidity). One can fit the satura¬ 
tion vapour pressure of water vapour over pure liquid water (mb) in the 
range -30 to +40C reasonably well by a law of the form P (T g ) = 

23.4 exp[0.076 (T -293)]. The surface mass density then becomes 
8 

P„ n (T ) " r (0.0177) exp[0.076(T - 293)] kg m" 3 (53) 

wo o S 

For 50% relative humidity and T = 288 , Eq. (53) gives 6.05 x 10 3 kg m 3 

8 

which is comparable to the midlatitude mean water vapour mass density 
given by the U.S. Standard Atmosphere, 1976. 

Using Eqs. (52 and 53) the effective heights z^ and tem¬ 
peratures can be calculated for the water vapour bands 1, 2 and 9. 

Letting r = 0.5 and T * 288°K we find, using the data of Table III-2, 

s 

the heights also given in Table 111-2. In all cases we have assumed strong 

lines. The heights imply temperatures via an atmospheric temperature - 

height profile. For T ■ 288°K we go directly to the U.S. Standard 

8 

Atmosphere. The temperatures T^ corresponding to the heights z^ are 

shown in Table 111-2. We will discuss below the temperature profile to 

be used when T j 4 288°K . 
s 
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c. Carbon Dioxide Bands (3, 4, 5, 7 and 8 ) 

The expressions for required equivalent width, Eq. (47), 
holds equally well for C0 2 as do the expressions for equivalent width 
Eqs. (48 and 49), provided z^ is relatively close to the surface, i.e., 
within about a scale height or two. Since CO^ is well mixed in the 
atmosphere up to stratospheric heights, it has the same scale height 
as the atmosphere as a whole. Thus for weak lines we have from Eq. (48): 


/ —z/H — z./H 

e dz - 1.86 v c ^ Sl He 


(54) 


where ra (v ) is the mass (volume) mixing ratio of C0„, p is the 

-3 

surface mass density of the atmosphere = 1.225 kg m and H is the 


scale height of the atmosphere. For strong lines equation (49) becomes 


W i * 2R i 


/" 

path 


D e 
C0 2 M o 


-2z/H 


1/2 


1/2 “ Z i /H 

= 2R ± (1.86 v cQ H/2) e 


1/2 " Z i /H 

= 1.93 R ± (v CQ H) ' e 


(55) 


In band 4 where the C0 2 absorption is very strong Eq. (55) with H = 8,500 m 
yields an effective radiating level which is on the order of 18 km. At 
18 km the scale height is appreciably different from its surface value. 

Hence it is appropriate to let H-*-H' where H' is a value more appropriate 
to the expected effective radiating level z^ . Here we choose H' - 7,300 m 
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which is appropriate to the density change between the surface and about 
18 km. 

In band 4 using H' ■ 7,300 m , T g - 288°K , v CQ - 333 ppm 
and R 3 from Table III-2 we find that z 3 - 18.2 km . Using the 1976 U.S. 
Standard Atmosphere this leads to an effective temperature T ®J217°K . 

Bands 3 and 5 contain the wings of the strong 15p CO^ ab¬ 
sorption band. These two bands have been included as separate bands be¬ 
cause it is here rather than in band 4 (which is already very strongly 
absorbed) that changes in CO^ concentration will be important so far as 
radiation balance is concerned. Although water vapour is also an impor¬ 
tant absorber in bands 3 and 5, we have neglected it here in the interest 
of simplicity. We believe the errors introduced to be minor ones since 
the scale height of water vapour is so much smaller than that of CO^* 

Most probably the errors introduced lead to an overestimate of the 
temperature increase for CO^ doubling. 

In bands 7 and 8, CC^ absorption is presently very weak 
and could well be neglected. However, as discussed in Sec. II1A of 
this report, future increases in atmospheric CO^ content should make a 
greater impact in bands 7 and 8 than in band 4. This is because band 7 
(at 10.4y) is presently a nearly transparent part of the atmospheric 
window. Similarly in band 8, (X^ absorption near 9.4y will add to the 
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ozone absorption present there. In bands 7 and 8 the CO^ absorption 
is presently weak and hence will grow linearly with increases in atmos¬ 
pheric CO 2 content; whereas in band 4 the CO^ absorption is strong and 
will grow only as the square root of the CO^ concentration. In addition 
band 4 is nearly completely opaque already. Using (54) for weak lines 
with z = 0 and the parameters of Table III-2, we find a mean transmis¬ 
sion of T = 0.95 for band 7. 


d. Ozone and COq Band (8) 

The ozone layer near 25 km altitude is a strong absorber 

at 9.6p. Since ozone is distributed in a layer and the stratosphere is 

held constant in our model, we treat band 8 as radiation from the surface 

at temperature T with a mean transmission given by 
s 


* - 1 - l(W C0 2 + H 0 > /A V 


(56) 


where W„ A and W are the equivalent widths of the 9.4p C0 o and 9.6p 0« 
absorption bands respectively. 


Assuming the ozone layer to be homogeneous with column 

density <cp> l = 3.68 x 10 ^ kg m ^ , R = 1731.8 cm ^(kg m , and 

°3 

mean pressure ^p} ** 25.5 mb have from Eq. (49). 


W Q - 2R q «cp>*) 1/2 «P>/P Q ) 1/2 


33 cm 


-1 


(57) 
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The equivalent width of the 9.4y CO 2 band is calculated as with the CO 2 
in band 8 yielding W « 7 cm ^ for v * 333 ppm , H - 8500 m , 

t»v2 CO 2 

-1 -2 -1 

and S C q * 1.2 cm (kg m ) at T ■ 288 K . Thus for the conditions 
of Table I1I-2 T ■ 0.6 for band 6. 


3. Temperature Profile 

The temperature profile T(z) is basically that of the U.S. 

Standard Atmosphere 1976. For a surface temperature T of 288°K the 

s 

temperature profile is thus the solid line in Fig. III-l. For values 

°f T g other than 288°K the lapse rate is maintained constant at 6.5 K/km 

as is the tropopause temperature (217°K). The structure above the tro- 

popause is assumed to remain that of Fig. Ill—1 regardless of surface 

temperature. Profiles for T + 288°K are illustrated by the dotted 

s 

lines of Fig. III-l. 


4. Comparison and Observations 

We now compare the IR fluxes predicted by our model with 
satellite measurements of the IR radiation emerging from several dif¬ 
ferent geographical areas. One set of such measurements - made by 
Hanel, et al. (1971) - from Nimbus 4, is shown in Fig. III-3. The 
thermal emission of the surface plus atmosphere emitted vertically 
upward was measured over a 95 km diameter circular area with a 
spectral resolution of 1.4 cm In each of the three graphs of 
Fig. III-3 the radiances of black bodies at various temperatures are 
shown for comparison. Also shown are the spectral bands of our model 
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Figure 111-3 

THERMAL EMISSION FROM THE EARTH'S SURFACE AND ATMOSPHERE 
EMITTED VERTICALLY UPWARDS AND MEASURED BY HANEL et al. 

(1971) USING THE INFRARED INTERFEROMETER SPECTROMETER ON NIMBUS 
4, (a) OVER SAHARA, (b) OVER MEDITERRANEAN AND (c) OVER ANTARCTIC 


The radiances of black bodies at various temperatures are superimposed as dashed lines. Average 
radiances predicted by the seven band radiation model discussed in the text are given by the bar graphs. 

In each case model radiances are given for three relative humidities. The solid-lines correspond to 

a relative humidity of 50% (r=0.5). The dashed-and dotted • • • • lines show changes for relative 

humidities of 10% and 90% respectively. The normalization factor a=1. 
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and the radiances predicted for each band. Model results are shown 
for relative humidities of 10, 50, and 90% (r - 0.1, 0.5 and 0.9). 

Here we have used v * 326 ppm appropriate to 1970 when the ob- 

CU2 

servations were made. 

For the Sahara desert case of Fig. III-3(a) we let T * 320°K 

s 

corresponding to the black body temperature in the atmospheric window 
(band 6). The effective heights and temperatures , determined 

as above, are shown in Table III-3. The nine band model gives a reason¬ 
ably good fit to the observational data. The agreement between model 
and observation is reassuring in that T g = 320°K is an extreme case— 
some 20°K higher than the highest zonal average shown in Fig. II1-4. 

For the Mediterranean Sea case of Fig. III-3(b) we again esti¬ 
mate the surface temperature from the atmospheric window (band 6) yield¬ 
ing T g = 283°K . A relative humidity of rather more than 50% seems 
appropriate as one could expect over water. Here again the model radi¬ 
ances are in reasonably good agreement with the observations. 

In Fig. III-3(c), results are shown for the Antarctic. Here 
the black body temperature in the atmospheric window indicates T g ~ 188°K 
(-85°C), very cold indeed. Since our temperature profile model calls 
for a constant tropopause temperature of 217°K, this indicates a 
temperature inversion and we assume a linear gradient of 5K/km from the 
surface to the 217°K level. The saturation water vapor pressure is 
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Table III-3 


EFFECTIVE HEIGHTS AND TEMPERATURES FOR THE SURFACE CONDITIONS OF FIG. II1-3 


(a) T - 320 K (b) T = 283 K (c) T - 188 K 

s s s 

Band Major r = 0.5 r * 0.5 r = 0.5 

i = Absorber z~ ± (km) T (K) z~ ± (km) T ± (K) z ± (km) T ± (K) 


1 

H 2 ° 

14.9 

223 

2 

h 2 o 

8.1 

268 

3 

co 2 

4.1 

293 

4 

co 2 

18.2 

217 

5 

C 0 2 

3.9 

294 

6 

None 

0.0 

320 

7 

co 2 

0 . 0 * 

320 * 

8 

° 3 /C0 2 

0 . 0 * 

320 * 

9 

h 2 ° 

1.2 

312 


10.3 

217 

0.0 

188 

3.4 

261 

0.0 

188 

2.0 

270 

0.0 

188 

18.2 

217 

19.7 

217 

2.1 

269 

0.0 

188 

0.0 

283 

0.0 

188 

0 . 0 * 

283 * 

0 . 0 * 

188 * 

0 . 0 * 

283 * 

0 . 0 * 

188 * 

0.0 

283 

0.0 

188 


Treated as optically thin. 
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Figure 111-4 

ATMOSPHERIC RADIATION FLUX TO SPACE F(T # ) 


Curves are given for the 9 band model for - 333 ppm and relative humidities of r * 0.1,0.5 and 0.9 
(solid lines). A model curve is also given for CO 2 doubling to Vqq^ = 666 ppm and r = 0.5. For comparison 

we include 3 curves (-and — lines) which represent straight line fits to observational data. The model 

curves are normalized to the average of the emperical curves at T $ •* 28Cf > K, i.e., a = 0.949 in Eq. (45). 
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so low at this low temperature that even high values of relative humidity 
still do not cause significant water vapor absorption. Thus the effec¬ 
tive radiating levels In bands 1, 2 and 9 are at the surface. The en¬ 
hanced radiance In the 15p CO,, band arises because the temperature at the 
effective height (In our model) Is actually greater than the very low 
surface temperature. Again the comparison shown in Fig. 1X1-3(c) is reassur¬ 
ing in that the agreement between model and observation is reasonably 
good even in this very extreme case. 

5. Radiated Flux 

a. Average Atmosphere 

As a first approximation, consider the Earth to have a 
homogeneous average atmosphere. This assumption allows us to extrapo¬ 
late conclusions about an atmosphere with a single temperature profile, 
etc. to the Earth as a whole. Clearly the climate model investigated 
in Sec. IV of this report will yield more precise results, but this 
initial approximation is still worth making. For our average atmos¬ 
phere we take the U.S. Standard Atmosphere, 1976 shown in Fig. 111-1. 

Using this profile and an average relative humidity of 50%, Eq. (46) 
yields the sum of the average radiated flux in each of the nine bands. 

The IK flux in each band as well as the effective heights and tempera- 

_2 

ture T^ are shown in Table III-2. The total IR flux F ■ 242 Wm when 

a * 1 . This value of F is slightly too high to balance the incoming 

-2 

solar flux F^(l-A)/4 • 235 Wm calculated when the albedo A has the 
conventional value of 0.33. However, more recent work indicates A is 
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-2 

about 0.29, yielding an incoming flux of 249 Wm . To make the out¬ 
going flux from our assumed average atmosphere balance the incoming 
2 

flux of 249 Wm we must let a * 1.03 . Having a > 1 Is to some ex¬ 
tent justified In that some compensation should be made tor Ifi flux which 
is not included within the 100-1230 cm * wavenumber band of the model, 
e.g., see the 1230-1500 cm * portion of Fig. 111-3. Another way to 
make the incoming and outgoing fluxes balance is to let a ■ 1 , but 

raise T from 288°K to about 290°K. 
s 

Although a is close to unity, it is not clear just what 

value a should take. From Fig. 111-4 we note that the data of Oerlemans 

and Van den Dool (1978) indicate F is about 238 Wm 2 at T • 288°K 

s 

based on zonally averaged fluxes and temperatures. These data also 

_2 

indicate a difference of about 4 Wm between the northern and southern 

-2 

hemispheres. If we take F as 238 Wm , a - 0.98 . 

b. Variations in F with Changes in and CO^ Content 

The key item of interest so far as climate effects are 

concerned is the alteration of F by changes in atmospheric constituents 

such as water vapor and CO^,. Consider the average atmosphere discussed 

above where v ■ 333 ppm , T g ■ 288°K , relative humidity • 50%, 

2 -2 

a * 1.03 and hence F “ 249 Wm balancing the incoming solar flux 
(Insolation). If the average relative humidity is reduced to 20Z apd 
other factors remain constant, F increases by 8.2 to about 257 Wm 2 . 

To regain balance with the incoming solar flux, the surface temperature 
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T would have to drop by about 3.9°K. If the average relative humidity 
s 

were to increase to 80%, F decreases by 0.9% and a surface temperature 
rise of about 1.1°K would be required to regain balance with the insola¬ 
tion. 


Again consider the aforementioned average atmosphere with 

= 288°K , r = 0.5 and a * 1.03 . Doubling the carbon dioxide con- 

-2 

centration from 333 to 666 ppm reduces F by 3.1% to 241 Wm . To rees¬ 
tablish radiation balance with the atmospheric CO^ content doubled, the 
surface temperature must rise by 3.8°K. A change in surface temperature 
could also imply a change in average relative humidity—probably an in¬ 
crease. If we increase the average relative humidity in our model from 
50 to 60% in addition to doubling the CO^ content, the new radiation 

balance is established at T ■ 290°K , an overall 4.3°K increase. If, 

s 

on the other hand, the relative humidity increases to 40% in addition 

to C0» doubling, the new radiation balance would occur at T * 288.2°K , 
z s 

an overall increase of only 2.8°K. Thus changes in relative humidity 
which occur concurrently with changes in CO^ content can serve to sig¬ 
nificantly enhance or diminish any CX >2 induced surface temperature 
changes. 


The surface temperature increase predicted by our model 
for a doubling of atmospheric CO^ content, namely 3.8°K, falls just 
above the range of 1.5 to 3°K suggested by Schneider (1975) as the 
range of reasonable values. He arrived at this range after reviewing 
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a number of model predictions. More recently he has suggested a range 
of 1.5 to 5°K which encompasses our result. 

We note that the change in outgoing infrared flux in¬ 
duced by a doubling of CO 2 in our model is caused by changes in bands 3, 
5, 7 and 8, i.e., the wings of the 15y band and the weak 9.4 and 10.4p 
absorption bands of CO^. Doubling the CO^ content causes the effective 
height in band 4, containing the core of the 15p CO^ band, to increase; 
but since the effective height remains in the tropopause region the 
effective temperature remains at 217°K and the flux radiated in band 4 
remains unchanged. 

c. Variations in F with Changes in Surface Temperature for 
Various Atmospheric Conditions 

By summing the radiation flux to space in each of the 

nine bands according to Eq. (45) we arrive at a total flux F . In 

Fig. III-4 we have plotted F as a function of surface temperature (T ) 

s 

for several sets of atmospheric parameters namely average relative 

humidity and carbon dioxide concentration. For comparison we show 

three emperical curves for F(T ) , which are straight lines fitted 

s 

to observations. We have normalized the model curves by letting the 
model curve for r = 0.5 and v m 333 ppm pass through the average 
of the emperical curves at T g * 288°K , the Earth's average surface 
temperature, i.e., a * 0.949 in Eq. (45). 
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Perhaps the most striking feature of the model curves 

themselves. Is the large effects caused by changes in average relative 

humidity (r) , especially for r £ 0.5 and T £ 270°K . This serves 

s 

to emphasize the importance of water vapor in the atmosphere and points 
out that changes in average relative humidity which might occur along 
with changes in CO 2 content could significantly alter predictions of (X^ 
climate effects based on an unchanging water vapor content. The atmos¬ 
pheric radiation model considered here includes the increase in satura¬ 
tion water vapor pressure with temperature; hence, for a constant rela¬ 
tive humidity, a CO^ induced increase in surface temperature causes 
an increase in water vapor content and hence an additional surface 
temperature increment over and above that caused by CO^ above. 

We note that for T <, 280°K , the three model curves for 

s 

F(T ) converge. This occurs because the decreasing temperature decreases 
s 

the water vapor content of the atmosphere for a given relative humidity. 

At low enough temperatures the water vapor absorption becomes negligible 
in bands 1, 2 and 9 and the model predicts the same flux regardless of r . 

An interesting feature of the model F(T ) curve for 
r * 0.5 and v * 666 ppm is that the difference relative to the 

CU2 

r * 0.5 , v * 333 ppm curve is most pronounced at higher temperatures. 

CU2 

This results from the temperature dependence of S in Eq. (54) for bands 
7 and 8. S increases by more than a factor of two between 250 and 300°K; 
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hence the greater CO^ absorption at high temperatures. This dependence 
of S on T in bands 7 and 8 would tend to mitigate CO^ induced tem¬ 
perature changes at high latitudes. 

Comparing the model curves with the observational curves 

we find that our model generally lies between the sets of observation 

data. In terms of the slope of the model curves, which is really the 

relevant comparison, the model curves for r £ 0.5 agree best with the 

Oerlemans and Van der Dool (1978) data at T £ 280°K and with the North 

s 

data for lower temperatures. The model curve for r = 0.1 is more closely 
in agreement with the Oerlemans and Van der Dool data over the entire tem¬ 
perature range. 

d. Linear Approximations to F(T ) 

In connection with energy balance climate models it is con¬ 
venient to approximate the atmospheric radiation flux to free space by a 
linear law of the form: 


F(T X ) « A + BT l (58) 

s s 

where T' is the surface temperature in degrees centigrade. This law 
s 

allows analytic solutions to a simple energy balance climate model as 
discussed elsewhere In this report. To provide an input to such models, 
we have used our nine band model to calculate A and B for varying CO^ 
content with relative humidity held constant at 50% and a * 0.949 . This 
value of a is the same as was used for normalization in Fig. III-4. 
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Table 111-4 gives the values of A and B for CC^ concentrations rang¬ 
ing from 300 to 1000 ppm. Comparing our values for A and B with 
those used by other authors, we find A - 209.9, 200.0 and 199.3 and 
B = 1.57, 1.45 and 2.40 for North (1975a), North (1975b) and Oerlemans 
and Van den Dool (1978) respectively. Comparing these values for A 
and B with our nine band model with v = 333 ppm and r ■ 0.5 we see 

from Table III-4 that our model values fall near the center of the obser¬ 
vational values. 

6. Summary and Conclusions 

Our hope was to include enough physics of the Earth's atmos¬ 
phere and of the radiation processes involved to construct a reasonably 
useful model. The extent to which this goal was achieved is shown graphi¬ 
cally in Figs. III-3 and III-4 where we compare model prediction with 
observational data. We conclude that our nine band radiation model, 
while far from complete, is useful in making at least crude estimates 
of the impact of increases in atmospheric C0^ content on future climate 
especially when coupled with the energy balance climate model discussed 
in section IV of this report. 

By assuming the Earth's atmosphere to be homogeneous with an 

average temperature of T ■ 288°K and an average relative humidity of 

s 

50%, we can estimate the change in surface temperature induced by a 
doubling of CO^ content. Our result is a 3.8°K increase which falls 
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Linear Approximation for K(T') Centered on A°C, F(T* in °C) * A+BT 


v co 2 <pp " ) 

A (Wnf 2 ) 

B(Wm~ 2 C -1 ) 

300 

203.2 

1.69 

333 

202.6 

1.66 

400 

201.1 

1.66 

500 

199.4 

1.61 

600 

197.7 

1.60 

666 

197.1 

1.57 

700 

196.6 

1.57 

800 

196.3 

1.54 

900 

194.4 

1.52 

1000 

194.3 

1.49 


Relative humidity m 50% and a « 0.949. 
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near the middle of the 1.5 - 5°K range suggested by Steven Schneider 
alter a review of model calculations. The CO^ doubling has its impact 
in the wings of the strong 15p CC^ absorption band and in the weak CO^ 
absorption bands near 10.A and 9.4p (bands 7 and 8 in our model). In 
band A which contains the core of the strong 15p CO^ absorption no 
changes occur for a CCL doubling. In this case although absoiption in¬ 
creases drive the effective radiating height upward, this height is in 
the constant temperature tropopause both before and after doubling; 
hence there is no change in radiated flux. Because we have neglected 

water vapor in bands 3 and 5, we expect that our figure of AT * 3.8°K 

a 

(for the average surface temperature increase upon CO^ doubling) is 

an overestimate. We have done a calculation in which the CO^ in bands 

3 and 5 was completely neglected. This resulted in a figure of AT - 1.4°K , 

s 

which we would regard as an underestimate. Since we expect the error in 

the larger figure to be smaller, we would correspondingly expect further 

refinements in our simple band model to lead to AT * 3°K . 

s 

The above result of AT q * 3.8°K surface temperature increase 
for a doubling of atmospheric CO^ 1 b in substantial agreement with the 
AT y * 2.8°K deduced in Section ili-A above. This is interesting be¬ 
cause these similar results were obtained by quite significantly dif¬ 
ferent approaches. Hence it lends credibility to the estimate of about 
3°K as the average surface temperature increase in response to a doubl¬ 
ing of atmospheric C0 7 content. 
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Certainly this nine band radiation model could be improved. 

For example in the handling of wavebands where there are two or more 
important absorbing constituents such as band 3 where both water vapor 
and CC >2 are important. However, one could easily question the wisdom 
of further model improvements in view of the reasonably good agreement 
with observations (Fig. I1I-3) and in the face of other uncertainties such 
as the response of water vapor to changes in surface temperature. 

It is no secret that water vapor is presently the dominant 
constituent in determining the atmospheric radiation balance (see Fig. 1II-2). 
Hence changes in atmospheric water vapor content which occur along with 
increases in CO^ content are clearly important and could be the dominant 
factor in the effect of CO^ on climate. Our model can quantify this 
importance in the sense that we can calculate the changes in radiation 
balance due to given concurrent changes in both H^O and CO^ content. For 
example, if we assume that upon a doubling of the average relative 
humidity also increases from 50 to 60%, our homogeneous atmosphere model 
implies a change in average surface temperature of 4.3°K, a nearly 15% 
higher temperature increase than with the relative humidity constant at 
50%. 


As shown by the dashed and solid lines in Fig. III-4, the differ¬ 
ence in F(T ) caused by a doubling of C0 o (with r held constant at 0.5) 
s L 

is not constant with variations in surface temperature. The difference 
increases with increasing T g . This results from the temperature depen¬ 
dence of the C0^ absorption coefficient in the weak 9.4 and 10.4|J bands. 
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While our nine band model is deficient in some respects, e.g., 
cloud variations are neglected, it does incorporate a modicum of at¬ 
mospheric physics and agrees reasonably well with experimental data. 
Hence it provides some insight into the climatic effects of increased 
atmospheric CO^ content. It also provides a radiative loss function 
dependent on average relative humidity and atmospheric CC^ content for 
use in the energy balance climate model discussed elsewhere. 
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Appendix A to Section III 
VIBRATIONAL BAND STRENGTHS 


The vibrational modes of a molecule interact with electromagnetic 
radiation because there is a significant charge associated with some of 
its vibrating constituents. If electronic excitations are ignored the 
molecule may be described as a collection of atomic masses Mj , such 
with an effective charge , interacting among themselves through 
complicated potentials. A gas of such molecules has a frequency depen¬ 
dent dielectric constant which satisfies the Kramers-Kronig relation; 


R a e(o)) 




lm e(o)*) a)* dm * 

,2 2 

a) - <jj 


(A-1) 


[This statement, equivalent to the analyticity of the function t(z) in 
the upper half plane together with the condition, from physical considera¬ 
tions, that e(u>) - 1 as u -* ® , follows Just from the fact that the 

molecule cannot respond to an external electric field before such a 
field is present.] 


The response of the molecule to a sufficiently high frequency ex¬ 
ternal electric field is independent of the interaction between the 
atom-like constituents and 


Lim t (u>) 

(tf+CO 



(A-2) 
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when n is the number density of molecules. When the high frequency limit 
is taken for the RHS of Eq. (A-l) we obtain the sum rule: 


-rf-f / 

J o 


Im £ (u) 1 ) duo 1 


(A-3) 


The integrated line intensity over a molecular vibration band is 
expressed in terms of £ as 


/ 

B 


S - 2 / *0 

nc 


(A-4) 


In a gas 


1Z % 

Im £* = Im [1 + (R£ £-1) + i Im£] 


~ \ Im e 


Then, because the integration range in Eq. (A-4) is less than that 
of Eq. (A-3), we have the inequality 


S < 2tt 


2\ " Q i_ 

V 


(A-5) 


For = x^x mass of a carbon atom and = y^e the inequality 


(A-5) is 


S < 2 . 10 


-16 


25 - 


(A-6) 
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Only very rarely is the bonding between atoms in a molecule so ionic 

that y^ ~ 1. In such cases large dipole moments may be present of 

order e(l-2A) - 10 ^ (cgs) * 10 D , as is found in F£,K , C£K , BeK . 

In most bonds D ~ 1 corresponding to y^ much less than 1. Then for 

-1 2 

H-C bonds (where would be of order 10 ), or C-C bonds y^ << 1. 

Values of S very much greater than a few xlO ^ cm would not be 
expected except in molecules where the number of relevant constituent 
atoms is so great that very many different vibrations in separate parts 
of the molecule have the same frequency. 
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IV 


CLIMATE MODELS 


The radiative models discussed in Section III neglect the redistribu¬ 
tion of the sun's energy by the movements of the atmosphere and oceans. 

The dynamics of the coupled ocean-atmosphere system remain an unsolved 
problem. In this section, we explore models in which the atmosphere 
and oceans are dealt with as a unit and the temperature of the earth 
is assumed to be a function only of latitude. 


A. Energy Budget Climate Models 

We wish to evaluate the effect of increased atmospheric concentra¬ 
tions of C0 2 on the climate, by examining simple climate models, the 
so called energy budget models. 

The most elementary energy budget model, which we choose to call 
zero'th order, simply asserts: 

7t4 Q 

4 (1 - Z) ( 1 ) 


where 
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o = Stefan Boltzmann constant 

T = effective mean radiating temperature 
e 

Q = solar constant 
a = mean albedo of the earth. 

From computations made at JASON 1978, we believe the relation between 
effective and surface temperature is given by 




I 


t 


T = (1 + 3t/4 T 
s e 




where x is a measure of the opacity of the opacity of the atmosphere 
whose dependence on CO^ concentration is described elsewhere in this 
report. For the current atmosphere, x , computed from first principles, 
was found to be 0.75. Putting in conventional values of the remaining 
constants: 

a = 5.6 x l(f 8 W/m 2 (°K) 4 
Q = 1360 W/m 2 
a = 0.3 

we get the quite acceptable figures = 257°K , and T g * 288°K = 15°C . 

As a predictive equation, Eq. (1) suffers from an array of defects. 
Changing, for example, the solar constant alters a in an unknown manner. 
A change in CO ^ concentrations affects the mean surface temperature, and 


t 






t 


« 


1 ' 
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hence also a 


Its effect on the biosphere, and hence again u , is, at the 


moment, a book of seven seals, needing intensive further investigation. 

A first attempt to overcome these difficulties is present in First 
Order Models . In these, surface temperature T is permitted to vary 
with latitude 0 , or one may imagine, if one prefers, that T is mean 
azimuthal temperature at latitude 6 . Heat transport is introduced, and 
the governing equation has the general form 

V + I * ^ (1-a) * q(l-a) 


where now 

V * transport mechanism 

I = outgoing radiative flux as function of T 
a * (mean) albedo of earth at latitude 0 as 
function of 0 and T. 

Various modellings of I, P, and a are possible. The simplest para- 
metrization of 1 - a is as follows: (It is convenient to introduce the 
variable x * sin 0 , in place of 0 .) 

jo. 7 if T > -10°C 
1 - a * Z(x) • ' 

10.4 if T ^ -10°C 
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where Z(x) is a correction for zenith angle, 0.7 is the co-albedo of 
land-sea, 0.4 that of ice-snow. In this way, ice albedo feedback is 
incorporated into the model. 


Better modellings would take account of land-sea mix at given lati¬ 
tude, have smoother transition of albedo from ice to non-ice, worry aboui 
cloudiness, etc., but the above is adequate to reveal the general features 
of these models. 


The radiative flux I may be modeled by permitting Eq. (1) to hold 
locally; thus 



__a_ 

1 +' .75x 


T 


4 


where r can now be adjusted for CO^ concentration. 

It is convenient, but not essential, to modify I by linearizing 
about Celsium zero, to get I » A + BT , with T now in degrees Celsius, 
appears, in practice, that A and B are then retrospectively adjusted to 
better fit current climate data, thus incorporating current cloudiness, 
for example. The effects of such a procedure should be obvious to the 
reader. 


At the same time, there are large, and rather serious discrepancies 
in the literature concerning the magnitude of B , which significantly 
affect the sensitivity of the model to changes in the solar constant. 
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We owe Che first consideration of first order energy budget models 
to Budyko and Sellers, independently. Subsequent work has been, in 
the main, variants of the schemes they contrived. 

Budyko's choice of D was: 

V - y(T - T) , 

T is the mean surface temperature, and y a phenomenological constant to 
be adjusted later. 

Sellers (and others) used diffusive transport 

p * - 3* [ D<1 * * 2) S] 

where D is, according to who you are reading, a constant, a function of 

dT 

x , a function of x and T , or even a function of x , T and — . 

Some writers have used a linear mix of Budyko and Sellers type 
transports. 

For the purposes of our immediate discussion, we will take 

v * - D h [ (1_x2) S] " div grad T * 

where D is a phenomenological constant to be determined later. 
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The first order energy budget equation becomes 


-D £[(l-* 2 > £] + A + BT - qZ(x) 


( 0.7 if T > -10°C 

( 2 ) 

( 0.4 if T < -10°C 


Equation (2) includes Eq. (1) from two points of view, so we may 
expect Eq. (2) to give reasonable mean surface temperatures. On the 
one hand as D gets large (more energetic diffusion), T gets more nearly 
constant and Eq. (2) becomes a version of Eq. (1). On the other hand, 
if Eq. (2) is integrated from -1 to +1 , the diffusive term vanishes, 
and what remains is a version of Eq. (1). 

Equation (2) is not a bad differential equation except for the 
driving term, which is a function of both the independent and dependent 
variable. It is usually solved as follows. Assume T is a decreasing 
function of x for 0 <_ x _< 1 , and even. Postulate x^ , 0 _< x^ _< 1 , the 
position of the northern ice line where T(x^) ■ -10 . Now the differential 
equation is an ordinary inhomogeneous one, requiring no boundary conditions, 
which can be solved by variation of parameters, or eigenfunction expansions. 

With the solution obtained, Q and D must be adjusted so that 
T(x^) = -10 , and monotonicity and evenness verified. 

If one supposes that for current climate, x^ * 0.95 , and q * ^ = 340 
D is determined uniquely. With D so determined, and q still fixed at 
340, one finds three positions altogether of the ice line meeting the 
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assumed requirements of evenness and monotonicity of T(x) , to wit: the 
normal, a heavily glaciated, and all ice. The heavily glaciated is 
not stable. 

With D fixed by the current climate, one may also investigate the 
effect of changes in S . A decrease of Q on the order of 2 to 3% is 
enough to give only an ice covered earth as solution, still on the assump¬ 
tion of evenness and monotonicity of T(x) . 

The results quoted above are mainly extracts from North 8 * 9 , though 
they are in general accord with the results of others 1 An outstand¬ 

ing exception occurs in Ref. 6, where a different determination of the 
constant B in the radiative mechanism I leads to strikingly distinct 
conclusions. 

There is one consideration on the methodology above which we feel 
is important to further consideration of such models. It concerns the 
role of the constant D , which is determined by matching current climate. 
It would be particularly desirable to give an a-priori estimate of D , or 
at least to present a persuasive argument that D is indeed constant over 
a reasonable range of the other parameters present. Otherwise, any 
predictions from the model of climatic response to change of the other 
parameters are highly suspect. 

This remark applies, mutus mutandi, to phenomenological constants in 
any model. 
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We have some further comments on the hidden or implicit assumptions 
in the procedure for solving Eq. (2). There is the assumption of even¬ 
ness of T(t) , (i.e., symmetry about the equator), as well as the assumption 
that T falls as one moves from equator to pole. There is nothing in 
the differential Eq. (2) requiring these assumptions, and as soon as 
one foregoes them, some interesting and provocative possibilities emerge. 

Let us for the moment assume that T rises in passing from 
equator to pole, but is an even function. One may, as before, solve the 
differential equation and hope the solution fits our assumptions. With 
some limited computer work, we found one such solution with the same D 
as for current climate, with the ice line at 30° latitude, but with a 
solar constant of 1680. This kind of solution, which we choose to call 
an anti-climate, seemed at first glance to be merely an aberrant, of no 
particular physical meaning or significance. Subsequently, it was 
discovered, however, that during the Permian era there was extensive 
glaciation in equatorial regions, but none in Antarctica, and anti¬ 
climates took on somewhat greater significance. Our geologic knowledge 
of Permian glaciation is, unfortunately, confounded by uncertainty as to 
the location of the continents. But aside from their possible relevance 
for ancient climates, we believe there is a variety of other anti-climates, 
though we must confess that our limited computer effort has not to date 
found any with the current solar constant and acceptable values of D . The 
computation to date does suggest, however, that there are many ice-band 
solutions, or ice-striped earths, in which one has a band of ice followed 
by a band of ice-free surface, etc. 
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There are, almost certainly, solutions of the differential equation 
not symmetric about the equator. These may find their physical realization 
in ice-age conditions, when the southern ice line remained in place, 
while the northern moved equatorward. 

It is conceivable that better understanding of First Order Models 
may be attained by trying to match well-documented ancient climates. 

We have not investigated stability of any of the bizarre climates 
alluded to above. We might add, however, that the possibility of bifurca¬ 
tion of the ice line does raise a question about the extant stability 
analysis of conventional climates. 

One of the most surprising features of First Order Models is how well 
they match current temperature profiles with such naive modelling of trans¬ 
port, radiation, and incoming energy. See, for example, the JASON Climate 
Model data elsewhere in this report. While this is a very pleasing 
development from one point of view, it is quite unfortunate from another. 
What it means, of course, is that only marginal evidence will be avail¬ 
able for the correctness of more sophisticated modelling of the important 
mechanisms. Bad further modelling may, nevertheless, reveal its nature 
very readily. 

We will address the question of testing more sophisticated effects 
in First Order Models presently. For the moment, we want to comment on 
the level of modelling incorporated in, say, Eq. (2). The I-R radiative 
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flux I should be nailed down as much as possible, and the constants 
A and B determined with greater precision. The modelling of the albedo 
can probably be improved; and we want to emphasize that, which is generally 
conceded by researchers, the effects of cloudiness must be assessed and 
incorporated into both the albedo and radiative flux. Since the relation 
between cloudiness and ground tempearture is unknown, speculations and 
investigations on the nature of cloudiness may, when incorporated into 
these models, produce feedback mechanisms with dramatic effects, or with 
very stabilizing effects. 

A word or two about the transport mechanism is in order. One of the 
more interesting features of these models is how little they appear to 
depend on the choice of transport mechanism. This is probably partly 
due to the freedom in choosing phenomenological constants, but if really 
so, should be explained on a-priori grounds. If the affirmative argument 
cannot be adduced, then efforts should be made to separate the distinct 
transport mechanisms on the basis of more sophisticated modelling, or 
any other reasonable tests. 

Finally, we want to indicate where we think these other tests may 
come from. Clearly, what is needed is data on other experiments which 
nature has performed for us. An obvious possibility, already mentioned, 
is the well-documented ancient climates. A richer source, probably, 
will come from examining seasonal and diurnal variations. A time 
dependent first order model will take the form 
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where C 


is an unknown heat capacity. An estimate of C may be obtained 
by simply looking at a time dependent zero f th order model 


C + a + BT * q(l - a) 

at a seasonal level, where T is mean temperature in the Northern hemis¬ 
phere, as a function of time of year, and q and a are known periodic 
functions of time of year. For each choice of C , the unique periodic 
solution may be compared with known seasonal temperature data in the Northern 
hemisphere to arrive at an optimal selection of C . 

The constant C so determined may then be utilized in a time dependent 
first order model: 


c + I + V =» q(l - a) 

to be analyzed first at the level of seasonal variation. The right hand 
driving terms would be parametrized as known functions of x , t , and 
the unknown periodic x^(t) . It is easy to Imagine iterative processes 
converging to the correct solution. One could, for example, postulate 
x^(t) as constant, solve the differential equation, use the solution to 
get another version of ice-line as function of time, and iterate. 
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In order to get realistic results, it may be necessary to build 
lapse times for melting of ice and freezing of water into this model. 

In general, it is quite difficult to see the consequences of pursuing 
this area of research, yet it is easy to see that it might be extraordi¬ 
narily fertile. 

As a next step in a research program, it would be useful to explore 
time dependent first order models at the level of diurnal variation. 

Beyond this lie second order models with a two-dimensional earth, 
and third order with an atmosphere. At the end of the line we have the 
G.C.M.s. 


B. JASON Climate Model 


To assess in a quantitative fashion the effects on climate of in¬ 
creasing the CO 2 concentration we have constructed and analyzed an energy 
balance climate model of the kind described in the previous section. 

The ingredients in an energy balance model are specifications of 

1) Incoming solar radiation and absorption by land and 
water: Insolation 

2) Outgoing radiation by the earth and by the atmosphere 

3) Transport of energy poleward by the ocean currents, by 
latent heat and by sensible heat. 
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More precisely we write 


C -r— * Insolation - Outward Radiation - Transport 
o t 

where C is some "heat capacity" which sets the time scale for the 
response of the temperature T to changes in any of the quantities on 
the right hand side of this equation. In principle, C can be a 
function of longitude, latitude and even time. 

Let us look at the three basic energy quantities listed above: 


1. Insolation 

This is the easiest and in many ways the most well established 
of the ingredients in our climate stew. The form of insolation is 


^Z(x)a(x,x ) (4) 

4 c 

where Q is the solar constant 

2 

Q a: 1380 watts/m 

and Z(x) is a zenith angle factor giving the mean annual distribution 
of insolation over latitude, x is the sine of the latitude and 


Z(x) - 1 - 0.482P 2 (x) 
= 1 - y P 2 (x) 


( 5 ) 
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1 2 

with P 2 (x) = (3x - 1) , the second order Legendre polynomial. The 

factor a(x,x c ) represents the absorption of the solar radiation by 

land sea and ice. We will use the idea of Budyko (1969) that a(x,x ) 

c 

varies with x for x < x c at which latitude the hemisphere is 

permanently ice covered and a(x,x ) x > x is a constant. Actual 

c c 

values for a(x,x c > are not precisely agreed upon by all workers, 
but we adopted the following values from North (1975b): 

x c = 0.95 , Latitude = 71.8° ; 


for x > x 

c 

a(x,x c ) = 0.38 , albedo - 0.62 (6) 

for x < x 

c 

a(x,x c ) = 0.697 - 0.0779P 2 (x) (7) 

which gives an average albedo of ~0.3 including land, water, and ice. 

Two comments are in order here: (1) we use zonally averaged 
quantities so that azimuthal or longitudinal dependence is washed out. 
This is an important simplification of the real two dimensional geometry 
of the earth's surface. However, following the orginators of energy 
balance models we are convinced it is an adequate starting point. 


4 > 


€ 


4 


4 
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4 
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(2) The particular values represented in Eqs. (5), (6) and (7) are, of 
course, important, but variations over small domains are unlikely to 
significantly affect our discussion or conclusions. 


2. Outgoing Radiation 

We desire a representation of the infrared radiation escaping 
the earth at each latitude as a function of the ground temperature . 

There appears to be agreement among workers that a linearized formula 

Radiation « + B^T(x) (8) 

is quite adequate. The coefficients and are, of course, inde¬ 

pendent of T(x) but may depend on cloudiness and other climatological 
features. Cess (1976) has reviewed the situation and his results 
derived from satellite observation are presented in Tables IV-1 and IV-2 

where A is the cloud cover fraction, 
c 

Table IV-1 

SUMMARY OF ZONAL ANNUAL CLIMATOLOGICAL DATA FOR THE NORTHERN HEMISPHERE 

The Units of F are W m - ^ 


Latitude 

(°N) 

T 

S 

(°C) 

A 

c 

F 

Observed 

Eq. (8) 

Error 

(%) 

5 

26.3 

0.51 

250 

252 

0.8 

15 

26.3 

0.44 

257 

259 

0.8 

25 

23.2 

0.41 

259 

256 

-1.2 

35 

15.9 

0.47 

241 

239 

-0.8 

45 

8.4 

0.57 

220 

218 

-0.9 

55 

2.2 

0.64 

204 

202 

-1.0 

65 

- 5.5 

0.64 

191 

190 

-0.5 

75 

-12.7 

0.61 

181 

182 

0.6 

85 

-18.0 

0.55 

179 

179 

0 

Average 

15.0 

0.51 

233 

234 

0.4 
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Table IV-2 

SUMMARY OF ZONAL ANNUAL CLIMATOLOGICAL DATA FOR THE 
SOUTHERN HEMISPHERE 
The Units of F are W m 


Latitude 

<°S) 

T 

(°6) 

A 

c 

F 

Observed 

Eq. (8) 

Error 

(%> 

5 

26.2 

0.50 

257 

263 

2.3 

15 

24.6 

0.47 

266 

263 

-1.1 

25 

21.4 

0.47 

262 

258 

-1.5 

35 

16.6 

0.54 

243 

245 

0.8 

45 

9.7 

0.65 

224 

225 

0.5 

55 

2.8 

0.79 

206 

202 

-1.9 

65 

-6.0 

0.77 

188 

190 

1.1 

75 

-33.0 

0.56 

163 

164 

0.6 

85 

-43.8 

0.47 

154 

154 

0 

Average 

13.4 

0.57 

236 

237 

0.4 


Using his values of average cloud cover we find for the Northern Hemis- 

2 2 

phere A = 210.59 watts/m and B = 1.57 watts/m -®C while for the 
R K 

2 2 

Southern Hemisphere ** 215.83 watts/m and B R - 1.59 watts/m -°C • 

In our actual calculations we adopted the values of North (1978) of 

2 2 
A= 211.2 watts/m and B R * 1.55 watts/m -°C . 


Now as discussed above, the variation of A and B with CO 

K K 

concentration or the concentration of any atmospheric constituent is 
given by 


ft 


\‘ o/ ( i+ t 'J 
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and 


B R- S/ ( 1 + ! T g) 


( 10 ) 


with T the net opacity at the ground and a and 3 some constants. By 
© 

using the relation between ground temperature, T , and the skin tempera¬ 
ture T needed to radiate back the net insolation 
e 


M-.) 


1/4 


(ID 


we determine that T s 3/4 for T = 257°K and T = 288°K - the global 

g e 

average temperature. 


As the opacity of the atmosphere to the infrared radiation 
entering in Eq. (8) increases, the radiation coefficients A R and B R will 
decrease and the average global temperature must increase. To get a 
feeling for this note that if we are in a steady state and integrate 
the energy balance equation over the globe it reads 

Infrared Radiation Out = Insolation 


As we change concentration, say, the net insolation is more or less 
constant (it is constant in our simple models if the ice line at 
doesn't change) so A R + B R (T) is constant. Is the global average 

temperature. Using the rules from Eqs. (9 and 10) we find, on the basis of 
the radiative transport theory presented before, that doubling the CO- 
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concentration from today’s value of 332 ppm would result in a change in 
global average temperature of 


■^Doubled C0 2 > - < T Today> ’ +2 ' 4 " C < 12 > 

This is a number consistent with a wide variety of other models, simple 
and complicated, and emerges because global energy balance really does 
represent the basic gross structure of our planetary temperature. Pre¬ 
cisely how this energy is transported about the globe does not affect in 
a dramatic fashion these global averages. Details of any climate model, 
including the one we are discussing here, will differ in how this 2°-3°C 
change is distributed over the globe. 



Following the ideas expressed by Stone (1973) with a somewhat 
different emphasis we have chosen to view this figure as representing 
several varieties of h eat transport cells . For example, suppose we 
consider the net heat flux R shown in Figure IV-1. It vanishes at 
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Figure IV-1 

THE AVERAGE ANNUAL LATITUDINAL DISTRIBUTION OF THE COMPONENTS 
OF THE ENERGY BALANCE OF THE EARTH-ATMOSPHERE SYSTEM IN KLY 
YEAR -1 (top) AND OF THE COMPONENTS OF THE POLEWARD ENERGY FLUX 

IN 10 19 KCAL YEAR" 1 (bottom) 
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x = 0 and x = 1 and peaks around x z 0.6 
by imagining a net heat flow which is 


We can parametrize this 



(13) 


where K is an "eddy diffusion constant" and the heat flux -Kx(l - x)dT/dx 

vanishes at x « 0 and x * 1 while peaking at x * 0.6 when dT/dx 

3/4 

behaves as * x . Really this is a representation of the actual heat 
flow by a lumped diffusion law with one overall constant K which sets 
the scale of the flow. 


If one wanted to make a finer kind of cellular structure, it is 
easy to imagine constructing several cells corresponding to latent heat 
alone or ocean currents alone, etc. For example, consider ocean currents. 
The sea surface temperature, call it T q , is different from the ground 
temperature T we have been considering heretofore. The heat flow due 
to ocean currents in the Northern Hemisphere is poleward from the 
equator to about x = 0.94 . One could try a representation of the 
ocean currents as 

* K o dt [ K(0 - 94 • X) “dll (14) 

with K q another "eddy diffusion constant" setting the scale for the 
kind of heat transport being considered—here ocean currents. 
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For latent heat in the Northern Hemisphere we would assign 
an effective temperature T^(x) to each latitude. Next note that the 
latent heat flow is negative for 0 x < 0.37 with a peak at x * 0.17 
and positive for 0.37 _< x 0.94 peaking at x s 0.59 . For this flux 
we would write 


-K 


1L 


d 

dx 1 


(0.94 - x)(x - 0.37)- 


d h 

dx 


+ k ,t / 

J 2L dx 


dT T 


x(0.37 - x)- 


dx. 


(15) 


with the "eddy diffusion constants" K > 0, i * 1,2 . A rule relating 

XL 

or coupling T o , and T would have to be constructed to complete 

this scheme. 


The general zone lying between a _< x _< 3 has the flux 

-K^[< 6 -*)(x-.)£] . (16) 

Here we will take one heat transport cell for each hemisphere 
and one diffusion constant K for each hemisphere. 


The climate model we have arrived at for the Northern Hemis¬ 
phere, say, is 

C If - K -k | x(1_x) |^| + \ + B R T(x>t) “ A Z < x > a < x > x c > • < 17 > 
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We will refer to this and its generalizations as JASON Climate Models 
(JCM) . 


Our first concern is with the allowed steady state versions of 


(1A): 


-K [*(1 - x)g] + \ + B r T(x) = J z(x)a(x,x c ) . (18) 


In this set z = 1 - 2x which turns (18) into 


^£(1 - z 2 )^] + o(o + l)T(z) = - z(x)a(z,z c ) - aJ = -S(z) ( 19 ) 


with 


o(c + 1) = -B r /K 


( 20 ) 


or 


2 1 V K A 


( 21 ) 


This is the Legendre equation. heat transport cell with flux given 
by (16) becomes Legendre's equation with z - [(a + $) - 2x]/(3 - a)|. 
The solution to (19) is 


T(z) = P q (z) f dw Q q (w)S(w) 

/_1 


+ V‘> f 


dwP (w)S(w) + CP (z) + nQ (z) , (22) 
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where P^(z) and (z) are the usual and second kind of Legendre 
function of Index a . £ and r) are Integration constants. The function 
Q 0 (z) is logarithmically singular at z ■ ± 1 as is P^(z) at z - -1 . We 
must choose £ and n so the temperature T(z) is finite for all z . Finite¬ 
ness at z ■ 1 requires n ■ 0 , while at z * -1 it is satisfied only when 

£ - - ff C °2 f dwS(w)P o (w) . (23) 

-1 


To derive these results we need 


2 sin air I *a 


P_(z) cos air - P a ( 


■ z) ] 


(24) 


for z along the real axis between -1 and 1. The result for T(z) is 
finally 


T(z) 


2 sin iro 


'.<■>/ 


dwS(w)P (-w) 
a 


+ P_(-z) (* dwS(w)P (w)! 
a j o' 

z 


(25) 


The solution for T(z) inside any heat transport cell has just this form 
with the o - - 1/2 + i B R /- 1 and the transport coefficient 

for that cell and z = (a + 0 - 2x)/(0 - a) in that cell. 
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• i 


It is important to note that we need require no conditions on 
T(z) or its derivatives except finiteness of T(z) throughout each cell. 
The "boundary conditions" discussed by North (1975a,b) and others are 
unnecessary. 


To determine K , the eddy diffusion coefficient in the JCM, we 
proceed by setting Q , A^ , B R and ■ 1 - 2x to their present day 
values given above. Then K is varied until T(z c ) achieves a definite 

2 

value, conventionally chosen to be -10°C. This yields a K * 1.37 watts/m -°C . 

2 

This diffusion coefficient is much larger than the value 0.59 watts/m °C 

found by North (1975b, 1978) when he considered a single heat transport 

cell running from pole to pole. Our heat flux proportional to x(l - x) 

is quite Inhibited in transporting heat near the equator, while North's 

2 

transport is proportional to J - x which is still sizeable near x ■ 0 . 

It is natural to expect the transport coefficient here to be larger 
since it is from the equatorial region that heat must be pumped. 


Once we have K such that T(z ) * -10°C we may calculate 

c 

T(z) for today's climate. This is given in Figure IV-2. 

This should be compared to the data recorded in Cess (1976) 
for both the Northern and Southern hemispheres (Figure IV-3). The JCM gives 
temperatures too low in midlatitudes, but rather reasonable at high and 
low latitudes. An effect like this should be expected from our dis¬ 
cussion of the structure of heat transport cells in the actual earth- 
atmosphere-ocean system. Because of the zeroes in latent heat transport 


« * 
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OBSERVATION DATA TAKEN FROM CESS (1976) FOR NORTHERN 
& SOUTHERN HEMISPHERES (0 is the latitude) 




and dips in sensible heat transport in raidlatitudes, less heat, in fact, 
is transported out of those latitudes than the simplest JCM would indi¬ 
cate. So a many celled JCM will push up the mid-latitude temperature as 
desired. 


The response of the JCM to changes in CO^ concentration is 

examined by changing A^ and according to (9) and (10). Using the 

results from earlier sections, we find a doubling of C0 2 in the atmo- 

2 

sphere to 664 ppm will lead to values of A,, = 207.7 watts/m and B = 

K. R 

2 

1.524 watts/m -°C . Since we have fixed our eddy diffusion parameter K , 

we adjust z so T(z ) = 10°C . We find that z = -1 for this 

c c c H 

and B , which means the ice line has retreated to the pole (see 
K 

Figure IV-4). So a doubling of CC >2 concentration according to the one heat 
transport cell/hemisphere JCM is likely to cause the ice caps to melt. 

The rate of increase of CC^ concentration has been 4.3%/yr 
of the industrial contribution. Taking 285 ppm as the preindustrial 
level we see that the fraction of CO^ in ppm increases as 

f__ = 285 ppm + (1.043) N • 47 ppm ; (26) 

t»u 2 

N is in years from today. So a doubling of C0 2 to 664 ppm should occur 
as soon as 50 years hence. If the ice line instantaneously follows the 
T = 10°C line as is assumed in energy-balance models, the ice line 
should retreat and disappear as C0 2 doubles. 
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Comparing Figures IV-2 and IV-4 we see that the average rise 
in T is slightly more than the 2.4°C evaluated from the conservation of 
AR + B R T . This is clearly because of the retreat of the ice line, and 
this is emphasized by the large rise in T in the polar region. 

These results are the essence of our analysis of the JCM. We 
have not explored to any extent questions such as the variation of the 
ice line as the solar constant varies or the sensitivity to the values 
of , B r , albedo, x^ , or the assumption that T(x c > * -10°C . These 
are interesting issues and important for one's confidence in the general 
scheme represented by the JCM. We feel it is rather clear that our heat 
transport cell description of heat flow in the earth-ocean-atmosphere 
system promises a realistic version of the actual complicated processes 
involved yet within the attractive simplicity of the energy balance 
models. 


An important issue we have not yet been able to explore in 
depth concerns the time dependence of the JCM. The problem is to study 
the solutions of Eq. (17) when , B R , x £ and possibly Q depend on 

time. This is important since the stability of the time independent or 
equilibrium states explored above is governed by the time dependent 
solutions. 


In (17) introduce 0(x,t) as 
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T (x, t) = exp 


cjf duB R (u) j 9<X - 


t) 


(27) 


so 0(x,t) satisfies 


I ~ K af [ x(1 " x) If] = [4 Z(x)a(x,x c (t» - A^t)] 

1 ^ 

x exp — I duB^(u) = S(x,t) . 


(28) 


# 




♦ 




Now go over to the Laplace transform of 0 


A 


T (x) 
P 



dte pt 0(x,t) 


(29) 


which satisfies 


-K 


_d 

dx 



Cpx (x) = S (x) - T(x,o) 
v P P 


(30) 


4 


4 


where 


S (x) = f dte pt S(x,t) 

P } o 


« 

(31) 


t 
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and we have noted T(x,o) * 0(x,o) . 


This equation is precisely of the form studied above: 



di (x)l 

x)—— I + n(n + I)t^(x) 


-S p (x) 


( 12 ) 


with 


n * 




P 


1 

4 


and 


S p (x) - [S p (x) -T(x,o)]/K . (33) 


So we know the solution, finite at z * ± 1 , to be (z - 1 - 2x) 


# 


* 


i (z) 
P 


2 sin irn 


I 


| P n <2 > 


C 


dwS (w)P (-w) 
P n 


+ P n ( - 2) / 


dwS (w)P (w) 

P n 


i 

i 


(34) 


♦ 


The temperature T(x,t) is recovered from 


T(x,t) - exp I 


if 1 1/ , o+I ” 

cf duB ii (u) ( f 

- 0 J' J a-i°° 


-i£ e PS (z) 
2irl e V ' 


) 


(35) 
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where Rea > 0 . This is the usual Laplace inversion formula and empha¬ 
sizes the point that the analytic structure of (x) in the p plane 
governs the stability of solutions to the JCM. Now what is attractive 
about this result is that given the temperature distribution at t ■ 0 , 
T(x,o) , the subsequent development of T(x,t) is governed by solutions 
to the same differential equation as the static JCM with the response 
dictated by the value of C . 

Furthermore the very important question of seasonal variations 
and their effect on accumulations of ice cover or cumulative effects on 
climate modification can now be explored. Seasonal variations are known 
for many of the ingredients to our model. North (1978) notes that the 
albedo factor entering z(x,x c > has a variation 

a(x,t) = 0.319 + P, (x)fo.071 cos 2irt/t 
1 *■ o 

+ 0.055 sin 2irt/t 1 + 0.202P_(x) , (36) 

o J 2 

with P^(x) = x and t Q * 1 year. Also the zenith angle factor entering 
the insolation behaves as 


Z(x,t) = 1 - 0.796 cop 2irt/t P. (x) - 0.477P o (x) . (37) 

o 1 2 


These ingredients can be straightforwardly folded into the time dependent 
JCM which may include as well a time dependence for the radiation para¬ 
meters and B R resulting from C0£ effects. The time scale of A^ 
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and B variations will be slow compared to t /2 tt a 1 month , so as a 

K O 

first approximation they could be taken constant. Detailed exploration 
of all this has yet to be carried out. 

One last note about T(x,t) in these models. The determina¬ 
tion of z c (t) must be done in a self-consistent fashion by looking at 

T(z c (t),t) - -10°C (38) 

Since T(x,t) is a functional of z (t) , this is a non-linear recursive 

c 

procedure which can in principle be carried out, though in practice it 
may be tedious. 


C. Sensitivity of Climate Models 

A normal procedure for assaying the degree of reality of a model 
in general and a climate model in particular is to test the sensitivity 
of the model to reasonable variations in its parametric makeup. 

In the following we use the notation of the preceding sections 
and for convenience we adopt the simplest energy balance model, i.e. 
that of North where the heat transport is given by 

- D fe £ «»> 
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Purely for consistency we keep this value of D constant, but again 
for convenience, following North we treat the radiated energy I as 
the dependent variable with a correspondingly different D which will 
be allowed to vary. 


In the expression for the absorbed solar heating 


Q[1+S 2 P 2 (x)] 


• b 

o 


X > X 

s 


• Ia o +a 2 P 2 (x)1 X < x 8 


( 40 ) 


we hold the parameters Q , b Q , a Q * a 2 

Since our problem is the sensitivity of climate to variation in CO^ 
of the atmosphere in the expression for the radiated energy 


A + BT 


(41) 


A and B will be varied in accordance with the discussion on the 
effect of CO^ on the radiant energy. This means that A and B will 
vary proportionately over a range of about 1% and D in the equation for 
1 will vary inversely. Increases in A and B will therefore imply 
a colder earth. 


152 


Digitized by 



We repeat the following initial values of the parameters: 


+ 



A 

211.2 

, 2 

watts/m 


B 

1.55 

watts/m^ deg C 

• 

D 

0.382 



b 

0 

0.38 



a 

o 

0.697 


# 

a 2 

-0.0799 



x„ 

0.95 



# The ice line temperature remains -10°C. There will be a slight 
inconsistency. North's value of D , 0.382, is only approximately, but 
closely, computed from the constants. Our computation will be exact 

# to the figures given. 

The approach is to perturb the constants A , B , and D and find 

# the new values of x , T (average T). and of course if needed, the new 

s 

temperature profile. What emerges is an extraordinary sensitivity of 
the ice line to slight variations in A and B . 
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The table is an indication of this sensitivity 


Table IV-3 


A/A 

O 

X 

s 

T 

1.02 

0.80 

5.45 

1.01 

0.87 

9.82 

1.00 

0.935 

13.81 

0.9975 

0.958 

14.92 

0.996 

0.980 

15.79 


A/A q of course is the fractional perturbation. 

The point of special interest is that the current climate corresponds 
to a sin latitude = 0.95, very close to ice free. A very small perturba¬ 
tion, 0.5% in the effect of CO 2 on the radiated flux, will drive the cli¬ 
mate to an ice free regime. 

In the other direction the area of the ice zone moves rapidly south¬ 
ward. 
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The curves shown in Fig. IV-5 are sample plots of the temperature 

variation for sample cases. For comparison an ice free regime A/A - 0.98 

o 

is shown although the formal solution for this case (x = 1) is trivial, 

s 

consisting as it does of the sum of terms with P q , P^ , anc * 

* 

P. of x as factors. 

4 

A number of obvious and curious conclusions can be drawn from 
these results. Perhaps the most interesting and possibly fruitful one 
is relative to the stability of the constants A and B . If the model 
has any validity at all this would mean that the real climate behavior 
must depend heavily on some specific processes that damp the effect of 
excursions in energy radiation or lock in the values of A and B in 
some fashion. 


This is why the 2-mode model works as well as it does. The assumed 
ice-line, y “ .95, represents a very small ice area: the exact solu^ 
tion for the ice-free situation is a good approximation to the "real 
case." 
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SAMPLE PLOTS OF TEMPERATURE VARIATION FOR SAMPLE CASES 






D. Sensitivity of Energy Budget Models to Changes in Albedo 


We briefly review the influence of such factors as albedo, infrared 
emissivity and cloud cover on climate. Although satellite data are now 
available there is some disagreement among various authors*^ ^ on how 
the data are to be interpreted. 


We first consider the arguments of Ref, 16, wherein ice-albedo feed¬ 
back is ignored. Using the global energy balance equation 

F - £ S(1 - a p ) (42) 


where F is the outward radiation flux, S the solar constant and ap 

the global albedo, we may write the "sensitivity parameter" 3 “ S dT /dS 

s 

in the form 


3 


F 


dF S 

dT 4 
8 



(43) 


where T g is the mean surface temperature. 


We note that in the case of constant albedo, 3 depends entirely on 
the IR feedback, giving 3 “ 70°C for an IR emissivity of .6, F * 233 W/m 
and T s = 288°K. The change of surface temperature due to a change in the 
solar constant of AS 


AT 

8 



(44) 
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is in this case . 7°C per 1% change in S (note for reference that 

dF/dT ■ 3.3) . Consideration of a more complete dynamics, with differing 
s 

feedback properties^produces 0 values from 35° to 400° C. 


The dependence of the infrared emissivity on surface temperature 
T g and cloud cover fraction has been parametrized by Budyko as 

F * 223 + 2.2 T - (47 + 1.6 T ) A (4:>) 

8 sc 


again with F in W/m^ and T in °C. 

19 

The most sensitive climates (0 s 400°C) were found by Budyko and 
20 

Sellers ; however. Ref. 18 claims that such sensitivity is greatly 
decreased If correct satellite data are used. All authors cited here 
use zonally averaged variables. 

The planetary albedo is decomposed as 

u*otA+a(l-A) (46) 

p c c s c 


where a is the clear sky albedo and a the cloud cover albedo, 
s c 

Nobody seems willing to predict the temperature dependence of A 

c 

Assuming that and u g are independent of T g , 0 is 





dF 

»A 

c 


t (a -a ) 
4 c 8 


dA 

_c 

dT 


( 47 ) 
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and argues that except near the poles 6 is sufficiently small (for 

subsequently deduced bounds on dA c /dT g ) the cloud cover correction is 

very small. However, the value obtained by Cess for 6 is small by 

virtue of cancellation between two large numbers. Using his fit 

(F * 257 + 1.6 - 91 A £ for the northern hemisphere), along with 

a * .44, a * -0.18 gives 6 * 91 - 88.4 = 2.6. Budyko's fit, Eq. (45), 
c s 

gives (for T * 15°C) 6 * 71 - 88.6 * -17.6, indicating greater sensi- 

s 

tivity. Clearly the cloud cover effect is very sensitive to the model 
and is therefore to be regarded at present as unreliable. It should be 
mentioned that the influence of cloudiness on albedo is completely ignored 
in Ref. 18. These authors disagree with Cess on the dependence of F on 
T g even though the same satellite data were used. The "correction" to 
Cess is held responsible for the climate being more stable to changes in 
S than claimed by Budyko and Sellers. 


The satellite fits cited above rely on comparison of different 

latitudinal zones. Thus the temperatures T are average annual values 

8 

suitable to different latitudes. Cess also discussed the seasonal tempera¬ 
ture dependence of F and a for a fixed zone, chosen to avoid compli¬ 
cations due to snow and ice (35°N). His motive for doing this was the 
apparent independence of the atmospheric feedback mechanisms on latitude 
(except possibly at high latitude). Again F and a are well fit to 

a linear dependence on T : 

8 
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f 


F = 185 + 3.5 T 

s 

a = 0.441 - 0.0086 T (49) 

p s 

These results are derived from monthly mean values of F and . If we 
use the values Eq. (49) blindly in Eq. (43) we get a very large value 
8 = 410°, very different from the value 8 = 145° obtained from Eq. (48) 
with negligible cloud feedback. 


f 


t 


♦ 


Cess proposes to interpret this difference as due to a strong 
positive seasonal feedback mechanism which does not occur for the 
long-term global climate change. His candidate for this is an effective 
temperature dependence of , to be estimated from the observed lati¬ 

tudinal dependence of a £ . Consideration of Eq. (46) then gives 


da 

_EL = 

dT 

s 



- A 

seasonal c 



(50) 


which restores agreement with 8 reckoned on an annual basis. 


* 


t 


A 


Finally, we mention Cess' estimates of dA /dT 

c 8 

find 



dT a -a dT 


s css 


From Eq. (46) we 

(51) 


Using Eq. (50) one finds values of -0.013, -0.020, -0.014 (°C)for 
25°N, 35°N and 45°N latitude, suggesting that the annual cloud cover 


A 


t 


A 


160 


Digitized by 


Google 



In this and other matters it is not 


decreases with increasing . 
clear whether latitude and temperature dependence have been correctly 
disentangled. 

Ice-albedo feedback was re-evaluated in Ref. 17 and was found 

19 20 

to be substantially smaller than originally thought. * 
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V RESEARCH NEEDS 

A CRITIQUE OF "A COMPREHENSIVE PLAN FOR 
CARBON DIOXIDE EFFECTS RESEARCH AND ASSESSMENT' 1 


A. General 


In May 1978, the Department of Energy's Office of Carbon Dioxide 
Effects Research and Assessment published a plan for research into the 
Global Climate Cycle and Climatic Effects of Increasing Carbon Dioxide. 

The Plan was the basis of a number of Research Program Development 
Papers prepared by specialists in various areas of carbon dioxide research. 

Our comments on the plan derive from our own research efforts as well as 
many discussions with scientists interested in the carbon dioxide question. 

The Department of Energy's plan is comprehensive and represents the 
results of a vigorous and wide-ranging effort to assess the opinions and 
ideas of those experts who have worked in the general area. It is useful 
as a compendium of possible research activities, but it can hardly be 
termed a research plan since the priority recommendations do not provide 
the overall guidance for the total effort. In particular, the plan does 
not indicate how the evolution of the program will depend on the results 
derived in any particular research effort. For example, until some rough 
estimate of the contribution of the biosphere is made, it is difficult to 
assess how much effort, if any, should be devoted to obtaining further detailed 
information on the biosphere. 
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The goal of the program should be to obtain a sufficiently good 
understanding of the carbon cycle and of the atmosphere to be able to 
predict with a high level of confidence climatic effects over time 
intervals of decades to centuries due to changes in atmospheric carbon 
dioxide concentration. Since there is not empirical evidence of the 
coupling between the observed Increase in atmospheric carbon dioxide 
and the current climate) we must turn to theory and indirect observations. 

With respect to an improved climate prediction capability, this 
difficult problem has attracted the attention of many outstanding 
scientists for a number of years, but progress has been slow. We 
believe that more effort along conventional lines will not lead to a 
timely resolution of the CO^ problem. The Plan does refer to simpli¬ 
fied climate models based on heat budget consideration, but does not 
give them a high priority. We believe that basically new approaches 
are needed, not necessarily based on simple energy budget models or 
complex General Circulation Models to deal with the CC^ question. 

The carbon cycle problem is more complicated. Some parts of the 
cycle have been studied for a long time, others only recognized 
recently. For example, the living carbon content of the oceans is 
relatively small, on the order of two Gtons, but the turnover time is 
short, about a week, far less than the terrestrial equivalent. The 
"box" models of the ocean tell us very little in that they are para¬ 
metrically insensitive. Much the same is true for the tracer approaches. 
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When a new one is Introduced, it carries with it its own set of 
uncertainties leaving the situation with one more equation, but also 
one more variable. The research program should recognize the limita¬ 
tions of eddy diffusive models and emphasize the content of the 
organic pools and turnover rates in both the ocean and terrestrial 
environments. 


B. Priorities 


Our set of priorities differs significantly from those given in 
the Plan. In our view, the primary problem in understanding the 
increase of carbon dioixde in the atmosphere is to determine the 
sign and magnitude of the biosphere contributions. If this contribution 
is close to zero as some chemical oceanographers would argue, then a 
good semi-quantitative description for changes in carbon dioxide is 
available and further research will only be a refinement. tf, on the 
other hand, the biosphere is as important as others have maintained, 
then we are far from understanding the sinks for the carbon based fuel 
injected into the atmosphere. Because of these considerations, first 
priority should be given to improved estimates of the contribution of the 
biosphere to recent changes in carbon dioxide, especially through the 
use of carbon isotopes. 


167 


Digitized by 



Our list of priorities, in order of importance is as follows: 




« 


• improved estimates of the contribution of the biosphere 
to recent changes in (X^, especially through the use of 
carbon isotopes; 


• detailed planning of oceanographic survey to assess the 
rate of transfer of the carbon species from the surface 

water to the deeper oceans; 4 


• investigations of more sophisticated energy budget 
models and statistical models of climate; 


• expansion of observing network to determine variation 
in the trace molecule composition of the atmosphere 
and the carbon isotope ratios; 


• feasibility studies of directly determining future 

changes in both vegetation and especially soil carbon; • 

and the ways these are being influenced by ground 
leasing and agricultural practices; 


t 

• investigation of the possibility of shallow-water 
sedimentation and dissolution contributing to the 
ocean CO^ balance; 

i 


• investigation of the seasonal variation of 00^ in the 
mixed layer in order to obtain better data on the 
turbulent diffusion parameters. 
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We would give an Improvement of the industrial source terms a much 

lower priority since these are known to within about 10%. We would also 

give research on the buffering factor a low priority. It does not make 

much difference at present just where in the range 1/7 to 1/11 the value 

lies. We do not believe that refinements of the gas-exchange rate across 

14 

the air-sea interface be given high priority. C measurements show that 
this rate is fast enough to insure spproximate atmosphere-mixed layer 
equilibrium on a time scale on order of magnitude less than that for CO^ 
growth or mixed layer residence times. The precise value does not enter 
into any prediction or interpretation of COj sources and sinks. 

In the area of improving climate prediction, the plan calls for 
support of coupled ocean-atmosphere General Circulation Models. We 
believe that the complexity of the models may be so great as to hide the 
underlying, physics and that the numerical methods employed may distort 
long term integrations. Research into models intermediate in complexity 
between the heat budget models and the GCM may provide useful insights 
into the prediction problem. 


C. Management Issues 

The research plan avoids questions of program management. We do 
have, however, serious reservations about housing program management 
within the Office of the Assistant Secretary of the Environment. The 
study of carbon dioxide is going to be a many-year effort of basic 
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research, involving a wide range of disciplines. The Office of the 
Assistant Secretary of the Environment is faced with a vast variety of 
day-to-day issues which may prevent adequate consideration of the longer 
term questions. The Office of Research has had a long history of managing 
multi-disciplinary long term efforts and is directly responsible for the 
National Laboratories. Because of their long involvement in basic research, 
the Office of Research would provide a more congenial home than the Office 
of the Assistant Secretary of the Environment. 

• Responsibility within the Department of Energy for the 
management of the carbon dioxide research and assess¬ 
ment program should be given to the Office of Research. 

The Plan refers to the international aspects of the carbon dioxide 
research program. We believe that the international nature of the pro¬ 
gram has received inadequate attention. Cooperation of other countries 
as well as such organizations as WMO and UNEP will be essential if the 
goal of the program is to be achieved. 

• International aspects of CO^ research should receive 
high priority. 
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